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SECTION  I 


INTRODUCTION 


The  fuel  distribution  achieved, by  transverse,  liquid 
injection  into  a  supersonic  airstream  (Fig.  1)  is  of  interest 
for  the  design  of  supersonic  combustion  ramjet  engines. 

In  prior  experiments  photographic  techniques  have  generally 
been  used  tq  measure  transverse  penetration  of  the  jet  bound¬ 
ary  for  both  subsonic  (Refs.  1-3)  and  supersonic  (Refs.  4-10)). 
gas:  streams,.  >'and/or  jet  boundary  lateral  spreading  (Refs. 

2,  5,  7,  and  9-11).  Ingebo's  (Ref.  12)  use  of  a  mass¬ 
sampling  probe  is  a  singular  case  in  which  in-stream  mea¬ 
surements  have  been  previously  applied  to  liquid  injection. 
These  past  studies  produced  several  empirical  formulas 
for  predicting  jet  penetration  (Refs.  1,.2,  7,  9,  and  .12) 
and  jet  spreading  (Refs.  2,  7,  9,  and  IT).  Theoretical 
studies,  however,  have  had  the  limited  goal  of  predicting 
jet  boundary  penetration  only.  Theoretical  models  genr 
erally  fall  into  two  categories:  those  that  assume  the  jet 
to  be  a  solid  body  (Ref.  6)  and  those  that  include  the  ef¬ 
fects  of  jet  disintegration  and  distortion  (Refs.  3,  8,  and 
13).  The  more  successful  models  depend  heavily  upon 
empiricism. 

The  existing  correlations  and  models  often  yield 
conflicting  predictions  of  the  absolute  and  relative  impor¬ 
tance  of  fluid  properties  such  as  injectant  viscosity  and: 
surface  tension.  In  fact,  the  ranges  of  conditions  tested 
in  many  experiments  are  insufficient  to  allow  the  evalua¬ 
tion  of  property  effects.  The  objectives  of  the  present 
program  were  to  advance  the  understanding  of  the  prob¬ 
lem  of  liquid  injection  by:  (a)  testing  through  wide  ranges 
of  system  properties,  (b)  using  more  comprehensive  ex¬ 
perimental  techniques  to  identify  and  evaluate  the  pro¬ 
cesses  and  properties  that  control  the  disintegration  of  a 
liquid  jet  column,  and  (c)  relating  these  processes  to  the 
resultant  spatial  distribution  of  the  injectant. 
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In  the  initial  phase  of  this  program,  emphasis  was 
put  on  the  determination  of  jet  boundary  location  by  photo¬ 
graphic  techniques.  Results  ^obtained  for  jet  penetration 
were  previously  reported  in  Ref.  14  and  are  only  briefly 
summarized  herein.  The  jet  spreading  results  are  pre¬ 
sented  herein,  but  they  are  also  discussed  in  Ref.  14.. 
Some  photographic  data  for  liquid  injection  obtained  with 
a  pulsed  ruby  laser  are  also  discussed  in  this  report. . 

The  more  recent  effort  irrihis  program  has  been 
directed  toward  s  determination  of  the  detailed  structure 
of  a- liquid  jet  flowfield  using  in-stream  instrumentation. 
In  particular,  a  mass-sampling  probe  was  fabricated  and 
evaluated,  and  together  with  pressure  probes,  applied  to, 
transverse  liquid  injection.  The  results. of  these  tests 
have:  (a)  led  to  better  understanding  of  liquid  dispersion 
in  a  supersonic  stream,  {b)  enabled  a  critical  evaulation 
of  existing  empirical  correlations  and  theoretical  treat¬ 
ments,  and  (c)  served  as  a  basis  for  an  improved  theo¬ 
retical  analysis,  as  described  herein. 
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SECTION  II 


EXPERIMENTAL  SETUP  AND  MEASUREMENT 
TECHNIQUES 

A.  TEST  MODEL,  INJECTANTS,  AND  TEST  CONDITIONS 

The  stainless-steel,  flat-plate  modeL(Fig,  2)  is 
0.  5 -inch  thick  and  5  inches  wide,  with  a  15°  leading-edge 
angle  and  a  10;  5-inch  working  length.  A  removable  injec¬ 
tor  insert  is  located  3.  06  inches  from  the  leading  edge  in 
a  plenum  that  is  fed  by  two  holes  drilled  180°  apart,  from 
opposite  sides- of  the  plate.  The  geometrically  similar  in¬ 
ternal  flow  passage  of  each  of  the  injector  inserts  cc  .sists 
of  a  cylindrical  exit  section  (diameter  =  length  =  dj)  that 
is  connected  by  a  transition  section  (blending  radius  = 

2  dj)  to  a  cylindrical  inlet  section  (diameter  =  5  dj).  Injec- 
tant  pressure  and  temperature  are  measured  in  the  plenum. 

The  liquid  supply  system  provides  for  remote, 
nitrogen-pressurization  up  to  1500  psi,  and  flow  rate  mea¬ 
surement  via  a  turbine  flow -meter.  Some  pertinent  prop¬ 
erties  of  the  liquids  are  listed  in  Table  1. 


Table  1 

Liquid  Properties  (at  20°C  unless  stated  otherwise) 


Liquid 

Density 

(gm/cm3) 

Viscosity 

(cp) 

Surface 

Tension 

(dynes/cm) 

Vapor 
Pressure 
(mm  Hg) 

Boiling 
Tempera¬ 
ture  (°C) 

Heat  of 
Vaporiza¬ 
tion  (cal/gm  ) 

Water 

0.  998 

1.005 

72.8 

37.54 

100.0 

539.6 

Isopropanol 

0.  785  > 

2.370 

21.  7 

-  22.07 

82.4 

159,4 

Refrigerant-11 

*  ;i 

1.464 

0.405 

19.  0+ 

672.  0 

23.8 

43.5 

at  boiling  temperature 
at  3 (PC 
*  at  25°C 


Preceding  page  Wank 


-  5  - 


For  water  and  isopropanol,  temperature -dependent 
properties  were  used  in  reported  calculations.  A  bother¬ 
some  feature  with  refrigerant- IT  is  that  the  injector  will 
become  -partially  plugged  withiceif  either  the  liquid:  sup¬ 
ply  system  or  the  air  supply  contains  water.  Extreme 
care  had  to:be  exercised::to?avoid  ice  formation,  which 
coulckbe  sensed^by  a,  variation  in  the. 'injectant  flow  rate 
.for  a  given  injectant  plenum  pressure. 

Table  2  lists  the  test  conditions  for  the  tests  dis¬ 
cussed  herein.  The  tests^  a  re  conducted  ina  2.  74-irich- 
,  ^diameter,  free-air  jet  from  a  contoured: supersonic  (Mach 
1.  62:  or  2j  72)  nozzle.  The  values  listed  for  the  effective 
injectant  orifice  diameter  de  are  based  on  calibration 
tests,  in  which  the  orifice  discharge  coefficient  was 
measured- for  Re -  between  4.  5  x  10?  and  1  x  l.ok,.  In  , those 
tests,  liquids  were  injected  into  a  quiescent  medium,  and 
the  metered  flow  rate  was  compared  with  the  ideal  value 
computed  by  the  Bernoulli  equation. 

B.  PHOTOGRA  PHIC  TECHNIQUES 

The  photographic  data  for  jet  penetration  were  ob¬ 
tained  by  direct-luminosity,  16-mm  cine-photographs  (31- 
ms  exposure).  Comparisons  with  data  from  70-mm 
schlieren  photographs  using  both*  20-ms  and  0.  25 -^s 
(spark)  exposures  showed  that  variances  in  measured  pene¬ 
tration  height  of  as  much<as,  40%  could  result  from  different 
photographic  techniques.  These  results  pointed  up  the  dif¬ 
ficulty  of  using  photographic  data  to  establish  any  weak  de¬ 
pendency  of  penetration  on  injectant  or  freestream  prop¬ 
erties  and  indicated  that  care  must  be  exercised  in  com¬ 
paring  results  from  different  experiments.  The  jet  spread¬ 
ing  data  were  obtained  by  direct-luminosity,  4-  by  5-inch 
still  photographs  (2.  5  ms  exposure). 

The  70-mm,  direct-luminosity  laser  photographs 
were  obtained  using  a  Q-switched  ruby  laser  having,  a 
nominal  rating  of  ~  150  MW,  with  a  pulse  width  of  ~  20  ns. 
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The  9/16-inch-diameter  beam  was  expanded  sufficiently 
to  give  complete  coverage  of -the  subject  and  the  correct 
exposure  for  the  Linagraph  Shellburst  (red  sensitive)  film; 
the  beam  was  not  recollimated.  To  decrease  the  effect  of 
ambient  light  without  the  use  of  extensive  shielding  pro¬ 
cedures,  the  laser  was  triggered  by  a  motorized  camera 
that  had  a  20-ms  shutter  window. 

C  MASS-SAMPLING  SYSTEM  AND  PRESSURE  PROBES 

The  main  requirements  for  a  sampling  system  for 
use  in  a  heterogeneous  stream  are  that  it  extract  a  repre-r 
sentative  sample  from  the  stream-  and  prevent  condensa¬ 
tion  and  entrapment  of  liquid  on  components  of  the  system 
during’  the  Collection  and  analysis  stages.  For  this  pur¬ 
pose,  -the  heated,  probe  shown  in  Rig.  3  was  fabricated 
and  evaluated.  Its  tip  consists  of  concentric  stainless 
steel  tubes  separated  by  an  aluminum  oxide  insulator.  The 
steel  tubing  serves  as  the  conductor  in  an  electric  circuit 
powered  by  a  12-volt  battery,  and  the  entire  length  of  the 
probe  and  tubing  to  the  sample  container  is  heated.  The 
probes  are  mounted' in  a  three-point  rake  with  a  0.  281- 
inch  center-to-center  spacing  (Fig.  4).  The  sample  con¬ 
tainers  are  wrapped  with  heating  tape,  so  that  subsequent 
analysis  can  be  made  by  gas  chromatography.  The  gas 
chromatograph  was  extensively  modified  to  assure  that  all 
internal  passages  are  kept  warm  enough  to  avoid  condensa¬ 
tion  within  the  unit. 

The  remaining  in-stream  instrumentation  consisted 
of  pitot-pressure  and  cone-static-pressure  rakes,  the 
most  recent  versions  of  which  art,  shown  in  Fig.  5.  The 
five-point  pitot  rake  used  0.  (165-inch  OD  by  0.  047-inch- 
IDsteel  tubing  spaced  at  0.281  inch  center-to-center.  The 
three  0.  J95-inch  OD,  cone-static-pressure  probes  are 
spaced  at  0.  562-inch  center-to-center.  Each  has  a  15? 
half-angle  cone,  with  four  0.015-inc.h-diameter,  equaCly- 
spaced  pressure  ports  located  0.  170-inch-from  the  tip 
and  connected  to  individual  strain-gauge  pressure  trans¬ 
ducers. 
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DIMENSIONS  ARE  IN  INCHES 


SILVER 

SOLDER 


304  STAINLESS  STEEL  TUBE 
0.109  OD  BY  0.012  WALL 


ELECTRICAL 

CONNECTIONS 


T^^gjii 

0.010 


DIAMETER 


ALUMINUM  OXIDE  TUBE 
0.078  OD  BY  0.008  WALL 


304  STAINLESS  STEEL  TUBE 
0.058  OD  BY  0.008  WALL 


'PROBE  RAKE  BODY 


FIGURE  3  HEATED  MASS-SAMPLING  PROBE  SCHEMATIC 
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FIGURE  4  TWO-PHASE  MASS-SAMPLING  PROBE 


PITOT  PROBE  CONE-STATIC  PROBE 


FIGURE  5  PHOTOGRAPH  OF  PRESSURE  PROBES 


SECTION:  III 

EXPERIMENT  A  L  RESULTS  FROM  PHOTOGRAPH  Y 


A.  SUMMARY  OF  PHOTOGRAPHIC  JET 'PENETRATION 
RESULTS  ?  v 

In  the  jet  penetration  studies  (ReL  14),  tests  were 
made  for  a  set  of  reference- conditions  that  were  arbitrarily 
chosen' as- follows : 


Water  injection:  d.  = 
3 

Airstream:  M 

00 


0.  0262  inch;  Tt  »-$7.(PR 
lj 


2.  72;  Tf  M  960°R; 

loo 


-p  =  14.  9  psia, 

‘CD 


A  series  of  tests  was  then  made  in  which  injectant  type, 

Ttj,  T^  dj,  and  M^  were  successively  varied;  in  most 

tests,  some  variation  of  pt  between «  300  and  «  900  psia 

j 

was  obtained,  _Our  results  (Fig.  6,  column  A)  showed  that 
the  parameter  yB/d^q0*  5  almost  collapses  the  data  for 
various  q's  in  a  given  test  to  a  single  curve,  and  the  de¬ 
viations  from  a  single  curve  have  a  random  nature.  More¬ 
over,  data  for  about  the  same  q  (Fig.  6,  column  B) 
showed  no  consistently  clear  effect  on  yB/de  of  any  other 
system  parameter  in  the  form  of  either  Reynolds  or 
Weber  numbers  or  physical  properties,  of  the  fluids.  When 
our  data  are  put  on  a  single  graph  as  in  Fig.  7,  the  data 
spread  is  nearly  constant  for  all  values  of  x/d e,  and  the 
following  equation  fits  the  data  to  within  ±  10%  in  the 
region  x/de  >  20: 

yB/de  =  1.  1  q°*  5  fn  (1  +  10  x/de)  .  (1) 

The  logarithmic  form  of  Eq.  (1)  was  suggested  jpy  two  cir¬ 
cumstances:  (a)  log-log  plots  of  yB/de  versus  x/de  showed 
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PENETRATION  PARAMETER 


DIMENSIONLESS  REDUCED  JET  PENETRATION 


T NEW  PREDICTION;  EQ;,(1) 

--  CHELKO  (1),EQ.  (2) 

- GEERYAND  MARGETTS  (2),  EO.  (3) 


- - —  CHELKO  (1) 

— -  DATA 


-  P—  q 

^  964  10.3 
4470  10.4 
4-591  10,3 


'! L  DIMENSIONLESS  AXIAL  DISTANCE  FROM  INJECTOR,  x/de 

FIGURE  7  COMPARISONS  OF  EMPIRICAL  PREDICTIONS  TO  PRESENT  DATA 


J  100r~ 


5 

8 

2  80 


CASTLE  (9),  EQ.  (5) 
EQ,  (6)  WITH  ft,  =  ptoo  (Moo  =  2,7) 

.EXPERIMENTAL, 

_  Pb  -  0'13  ptj,  EQ. 


5 

re 

l~ 

u 

z 

m  - 

CL 

2 

i  20 


^  EQ.  (6)  WITH  Pb  =  Poo 
■  KOLPIN,  ET  AL  (7),  EQ.  (4)  1 


JET-TO-FREESTREAM  PRESSURE  RATIO,  p,./poo 


FIGURE  8  COMPARISONS  OF  EMPIRICAL  PREDICTIONS  FOR 
PENETRATION  AT  x  =  150  TO  PRESENT  DATA 


that  significant  curvature  remained. in-the  trajectory  curves, 
and  (b)  the  logarithmic  form  often  results  in  theoretical 
studies  of  liquid  injection  (see  Section  V).  The  equation  is 
compared  to  the  jet  penetration  results  obtained  by  in- 
stream  mass -sampling  in  Section  IV. 


Our  jet  penetration  data  are  in  fair  agreement  with 
some  previously  derived  empirical  correlations.  The  fol¬ 
lowing  two  from  Refs.  1  and  2,  respectively,  predict  that 
yjj  varies  with  x  and  q  according  to  simple  povter  laws: 


.  ,c  0.  265  0.475  0.22 
0.45/9  q  x 


(2) 


0.5  0.27 


yB  =  2,0q--x 


(3) 


[Note  that  Eqs.  (2)  and  (3)  are  scaled  in  terms  of  d., 
rather  than  de;  hence,  Cd  =  1  is  assumed  in  the  compari¬ 
sons.  ]  Equations  (2)  and  (3  rare  compared  with  the  data 
imFig.  7,  where  the  two  curves  shown  for  Eq.  (2)  are 
based  on  the  extreme  values  of  p  and  q  that  existed  in  our 
tests.  In  view  of  typical  photographic  data  bias,  the  quan¬ 
titative  predictions  of  Eqs.  (2)  and  (3)  compare  favorably 
to  the  data.  However,  the  power-law  dependency  of  y-g 
on  x  gives  poorer  overall  qualitative  agreement.  In  addi¬ 
tion,  the  p  dependency  predicted  in  Eq.  (2)  is  not  supported 
by  the  present  data,  as  seen  in  the  lower  part  of  Fig.  7, 
which  compares  trajectories  for  which  p  was  varied  but 
the  q's  were  nearly  equals 


The  following  correlations  from  Refs.  7  and  9, 
respectively,  predict  that  maximum  yg/de  always  occurs 
by  Sc/tig  =  159  and  is  given.by: 


3 

(y„/d  )  =  (15,  7/M  )  (p,  /pj 

erniax  “  t. 

0 


0.51 


0.407 


(4) 

(5) - 
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In  order  to  compare  Eqs;  (4)  and  (5)  to  our  data,  which 
•  were  limited  to  x/de  <  100,  Eq.  (1)  was  evaluated  for 
x/de  =  150,  and,  using  the  Bernoulli  equation  and  the  de¬ 
finition  of  q,  rewritten  as  follows: 

(ys'Vmax  -  ^  62/M->  [<Pt.  +  pb 5  •  <6> 

3 

Equation  (6)  is  compared  with  Eqs.  (4)  and  (5)  in  Fig.  8, 
where  two  cases  of  the  back  pressure  sensed  by  the  jet 
are  considered:  p^  =  p*,  and  -  pt  ,  where  pj  is  the 
freestream  pitot  pressure.  These  t^o  cases  bracket  the 
experimental  pb's,  which  roughly  agreed  with  the  express 
sion: 

p^  «  0, 13  p^  .  (7) 

3 

This  latter  expression  allows  Eq.  (6)  to  be  written  as: 

(yB/de»max“<10'2/M.)(pt./p»)0'5  ' 

l 

Equation  (6a),  which  approximately  represents  the  data, 
also  is  shown  in  Fig.  8  to  fall  between  Eqs.  (4)  and  (5), 

It  is  in  closer  qualitative  agreement  (slopes)  with  the 
former,  but  is  in  closer  quantitative  agreement  with  the 
latter.  In  fact,  the  predictions  of  Eq.  (5)  are  50  to  70% 
higher  than  those  of  Eq.  (4)  for  the  range  of  conditions 
covered  in  our  tests.  It  is  apparent  from  the  compari¬ 
sons  of  Eq.  (1)  [which  lead  to  the  approximation  (6a))  with 

Eqs.  (4)  and  (5)  that  the  use  of  p^./p^  as  an  independent 

3 

variable  is  essentially  equivalent  to  the  use  of  q  for  the 
ranges  of  conditions  normally  encountered  in  liquid  injec¬ 
tion.  However,  the  frequent  and  natural  appearance  of  q 
in  analytical  studies  suggests  that  the  use  of  q  in  em¬ 
pirical  correlations  is  more  justifiable. 

From  our  data,  it  was  established  that  a  superior 
correlation  was  obtained  with  dg,  as  opposed  to  dj,  as  the 
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length,  scale.  This  is  apparently  due  to  the  fact  that  the 
use  of  de  partially  accounts  .for  the  boundary  layer  of  the 
emerging  jet,  whv'h  depends  upon  the  internal  geometry 
of  the  injector.  In  fact,  an  effective  diameter  based  on 
momentum  flow  considerations  would  be  even  more  appro¬ 
priate  than  the  present  de,  which  is  based  on  mass  flow 
considerations. 

B.  PHOTOGRAPHIC  STUDIES  OF  JET  BOUNDARY 
SPREADING 


Effects  of  d.:  (0.  017  arid  0.  027  inch),  p^.  (118  to  926 

psia),  and  M»(l.  62  and  2.  72)  on  the  spreading  of  isopro- 
panol  jets  were  determined,.  Still  photographs  of  2.  5-ms 
exposure  (Fig.  9)  were  obtained  by  recording  the  back- 
scattered  light,  from  flood  lamps  located  above  the  flat- 
plate  surface*  The  data  were  subject  to  variations  up  to 
25%„as  a  result  of  individual  interpretation  of  the  jet 
boundary  location.  It.  was  possible  to  make  measurements 
at  (zg/de)‘s  up  to  «  70  for  d.  =  0,027  inch  and  up  to^  125 
for  d-  =  0.  017  inch.  J 

V 


The  jtypical  results  in  Fig.  10  show  that,  for  given 
dj  and  M«,  Zg/ d'e  increases  with,  p^.  for  x/de  ^  10;  however, 

in  some  cases,  the  opposite  trend  is  observed  for  small 
(x/de)'s  (this  behavior  at  Mw  =  2.  72  can  be  seen  in  Fig.  10, 
where  the  leading  edge  of  ;che  jet  at  pt  =  118=psia  is  notice¬ 
ably  blunter  for  x/de  s  5),  For  given^d.  and  p^.  ,  zR/d 

decreases  with  Mm  for  x/de  >  3  but  increases  with  M,,,  for 
smaller  x/de's.  This  behavior  suggests  that  initial  jet 
spreading  i3  controlled>by  dynamic  and/or  viscous  forces 
that\the  airstream  imposes  within  the  stagnation  region 
of  the  jet  column,  tearing  ligaments  or  drops  from  the  jet 
and  imparting  lateral  momentum  to  them.  These  initial 
spreading  effects  are  greater  when  Mw  increases  and  when 
the  mean  vertical  momentum  of  the  jet  is  smaller  (lower 
p^.  for  given  d.)  or  when  the  jet  has  a  lower  Re.  and;  hence 
a  ^  relatively  larger  boundary  layer  (smaller  d^  for  given 
p^.).  Having  been  torn  away,  the  drops  are  subjected  to 
3 
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FIGURE  9 


(a)  P,  =  118  PSIA,  Moo  =  2.72 
i 


(b)  Pt.  =  925  PSIA,  Moo  »  2.72 


(c)  P,  =926  PSIA, Moo  =  1.62 
i 


0  20  40  60 

x/dj 


PHOTOGRAPHS  OF  JET  SPREADING  FROM  0.027-INCH  ORIFICE 
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DIMENSIONLESS  JET  WIDTH,  zB/de 


■H 

(a) ^EFFECTS  QFMNJECTION  PRESSURErAND 
FREESTREAM  MACH  NUMBER 


(b)  EFFECTS  OF  INJECTION  PRESSURE  AND  ORIFICE  DIAMETER 
FIGURE  10  TYPICAL  JET  SPREADING  RESULTS. 


-  20 


dynamic  forces  imposed  by  the  freestream.  Since  the 
resulting  axial  momenturh  interchange- also  increases 
with  the  ultimate  spreading  decreases  for  an  in¬ 
crease  in  Mw. 

The  only  substantial  efforts  to  correlate  jet  .spread¬ 
ing  data  have  been  directed  toward  predicting  the  maximum 

jet  width  Zt>  /d  ,  which  is  assumed  to  occur  by  x/d„  = 
j  -tJniax  e  '  J  e 

150.  (Note  that  the  change  in  character  of  jet  width  versus 
system  parameters  with  increasing  x  requires  a  rather 
complicated  expression  to  predict  the  entire  boundary  tra¬ 
jectory.  )  Initially,  the  correlations  were  in  the  form 
zBmax^de  *  (Ptj/Po,)0'  where  c  was  determined  as  0. 133 

and  0. 197  in  Refs.  7  and-9,  respectively.  Photographic 
data  were  obtained  in  Ref.  7  for  water  injection  into  a 
Mach  2,.  8  airstream,  with  de’s  of  0.  0117  to  0.  0211  inch. 
The  correlation  from  Ref.  9  was  based  on  photographic 
data  for  the  injection  of  water,,  heptane,  and  a  58%  glyc- 
erin-42%  water  solution  into  Mach  2.  8  and  Mach  4.  0  air 
streams  from  orifices  with  de!>s  of  0.  0117  to  0.  0624  inch. 

It.  Ref.  11,  the  correlation  was  modified  to  include  a  Mach 
number  dependency: 

ZB  M«  /de  =  45  (pt  •  (8) 

max  j 

The  present  data  are  compared  to  Eq.  (8)  in  Fig.  11, 
which  includes  some  typical  data  points  from  Ref.  9  for 
injection  of  water,  heptane,  and  a  water-glycerin  solution. 
To  make  the  comparison,  it  was  necessary  to  extrapolate 
some  of  the  spreading  data  curves  to  determine  zBmax^de’ 

However,  such  extrapolation  was  reasonably  safe,  since 
the  slopes  of  the  curves  in  Fig.  10  had  become  very  small 
at  x/de  >  60.  The  present  data  (open  symbols  in  Fig.  11.) 
are  represented  well  (±  7.  5%)  by  a  single  sti’aight  line, 
but  the  equation  for  this  line  is: 

ZB  M.0-5/de  =  10.4(pt/p„)°-25.  (9) 

max  j 
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MAXIMUM  LATERAL  SPREAD 
PARAMETER,  (fBMoo°-5/de)  m#! 


FIGURE  It  COMPARISON  OF  JET  SPREADING  DATA  TO  EXISTING  DATA  AND 
CORRELATION 


Equation  (8)  predicts  jet  widths  that  are  higher  than  the 

present  data  by  factors  of  «  2’ (at  =  1000),  to  «  3.  3 

J 


(Pt./p  =10).  However,  owing  to  the  considerable  data, 
3 

scatter,,  a  line  having- a  0.  25  slope  does  not  violate  the 
Ref,  9  data. 


Mass-sampling  data  at  x/de  =  51  for  refrigerant- 11 
with  M'ob  =  2.  72  that  are  discussed  later  yield  Zg  M#^  ^/de 

36.  3  for  pt./pw  =  61, 4,  ai?-»  zB  Ma°f  D/de  =  27.  9  for  pt./pa 
3  3 

20.  7.  From  the  M*  =  2.  72  results  shown  in  Fig.  10,  it  is 
judged- that  would  increase  by  no  more  than  10%  in  go¬ 
ing  from  x/de  =  51  to  the  point  where  ZBmax/de  occurs. 

Even  if  a  15%  increase  were  accepted,  upper  limits  for 
these  two  data  points,  which  are  plotted  in  Fig.  11,  would 

be  41.  7  and  32. 1  at  pt./pw=  61.4  and  20.  7,  respectively. 

3 

A  line  through  these  . points  still  has  a  slope  of  0.  25,  but  a 
constant  of  15.  0;  compared  to  the  10.4  in  Eq.  (9).  Rela¬ 
tive  to  this  dashed  line,  predictions  from  Eq.  (8)  would 
still  be  >  84%  high  for  the  range  of  the  present  data,  and 
as  much  as  67%  high  for  the  data  region. (p^/p^  >  150)  on 
which  Eq.  (8)-- is  based..  J 

If  the  procedure  used  to  derive  Eq.  (6a)  from  Eq. 
(1)  were  applied  inversely  to  Eq.  (9),  there -would  result: 
zB^de  w  25*  In  fact,  the  present  photographic  data  are 
equally  well  correlated1  (±  7.  5%)  by 

z_/d  “  10r  q°'  (10) 

AJ  V 

In  any  event,  the  correlations  found  for  yB  /d  and 

v#  2 

Zn  / d  combine  to  give  y-™  / d  «  0.  08  (z-r  /d  )  , 

Bmax  e  -  b  ;BmaX'  e  -“max  e 


C.  LASER  PHOTOGRAPHS  OF  LIQUID  INJECTION 

the  object  of  the  laser-photography  studies  was  to 
determine  whether  the  flowfield  could  be  sufficiently 
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stopped  to  allow  identification  of  individual  fluid  particles. 
Figures.'-' T2  shows  direct-luminosity  photographs  obtained 
with  the  camera  set  at  f /4.  5,  with  a  1.  06:1  magnification. 
Parts  a  and  b  show  results  for  identical  test  conditions  at 
different  times,  and  parts  c  and  d  show  results  for  a  lower 
arid  a  higher  injection  pressure,  respectively.  (The  probes 
seen1  in  the  photographs  are  not. in  the  jet  flowfield.  The 
reflection  of  the  jet  flow  is  also  seen  in  the  flat-plate  sur¬ 
face.  )5  These  pictures  3how  the  unsteadiness  of  the  flow 
and  the  irregular  shape  of  the  outer  jet  boundary;  These 
features  had^been  noted  previously  using  a  250-ns  spark 
source,  as  compared  to  ~  20  ns  for  , the  laser;  however, 
the  laser  photographs  show  considerable  more  detail 
within  the  jet  core,  arid  one  can  clearly  see  variations  of 
the  injectant  density  over  the  flowfield.  The  jet  remains 
essentially  intact  for  some  distance  before  its  breakup 
begins,  which  results  in  a  region  behind  the  jet  and  near 
the  plate  surface  that  is  nearly  devoid  of  injectant  (farther 
downstream,  fluid  mixing  eventually  eliminates  this 
region).  Unfortunately,  the  number  of  fluid  particles  is 
too  large  to  permit  therii  to  be  identified  individually.  It 
seems  doubtful  that  the  more  sophisticated  holographic 
techniques,  would  solve  this  problem  either. 
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SECTION  IV 

IN-STREAM  FLOWFIELD  MEASUREMENTS 

A.  EVALUATION  OF  MASS-SAMPLING  SYSTEM 

It  was,  recognized  early  in  this  program  that: 

1.  The  spatial  distribution  of  the  injected  mass 
must  be  known  if  one  is  to  characterize  completely  the 
injection  process; 

2.  Standard' photographic  techniques  are  not  sensi¬ 
tive  enough  to  permit  an  evaluation  of  the  effects  of  injec- 
tant  properties  such  as  viscosity  and  surface  tension  on 
this=  spatiai:  distribut ion;  a nd • 

3.  Measurements  of  this  spatial  distribution 
would  improve  the  understanding  of  the  controlling  mechanr 
isms  and  form  the  basis  for  realistic  theoretical  models 

of  liquid  injection. 

In  some  early  tests  a  refrigerant- 2 2  (chlorodi- 
fluoromethane)  jet  was  sampled  with  an  unheated  probe, 
and  it  was  found  that,  on  the  jet  centerline,  the  maximum 
injectant  composition  occurred  at  a  height  that  was  only 
40%  of  the  maximum  penetration.  The  use  of, an  unheated 
probe  was  permissible  in  those  tests  because  the  high 
vapor  pressure  of  refrigerant- 2 2  assured  that  it  would 
rapidly  vaporize  in  the  airstream  and  not  condense  in  the 
sampling  system.  To  test  less  volatile  liquids,  it  was 
necessary  to  develop  the  heated  sampling  system.  The 
probe  design  and  use  of  an  insulated  and  warmed  gas 
chromatograph  were  described  in  Section  IIC. 

To  evaluate  the  sampling  system,  premixed  sam¬ 
ples  of  known  composition  were  collected  and  analyzed. 

The  test  apparatus  (Fig.  13)  for  preparing  the  sample 


-  28  - 


employed  contra-stream  injection  of  a  metered,  liquid  flow 
via  a  spray  nozzle  into  a  metered  airstream  that  could  be 
heated  to  500°F.  Vaporization  and  mixing  occurred  in  a 
1.06-inch-diameter,  81 -inch-long  pipe.  The  mixture  is 
discharged  either  directly  from  the  pipe  as  a  subsonic 
stream,  or  from  a  1.  3-inchrexit-diameter,  Mach  2  noz¬ 
zle.  In  the  initial  tests,  subsonic,  isopropariol-air  mix¬ 
tures  .were  sampled  and  analyzed, as  being  more  than  30% 
lower  ,than  the  input  values.  These' results  led  to  a  deci¬ 
sion  to  use  the  more  volatile  refrigerant- 11  (trichloro- 
fluoromethane),  which  boils  at  75°F  and  reduces  heating 
requirements  on  the  sampling,  equipment. 

The  results  with  refrigerant- 11  are  summarized  in 
Table  3.  Tests  A-I  were  conducted  with  heated  air  to  avoid 
condensation  Of  the  mixture  on  pipe  wall  surfaces.  Except 
for  the  supersonic  Test  E,  the  sampled  mass  percent  was 
lower  than  the  input  value.  The  errors  for  the  subsonic 
tests  varied  from  -8  to  -17%  for  the  supersonic  tests, 
from  6  to  7%.  However,  samples  taken  from  different 
points  of  the  stream  usually  were  in  close  agreement, 
suggesting  thorough  mixing  of  the  stream.  No  consistent 
dependence  of  sampling  error  on  the  input  mass  percent 
level  was  observed. 

Table  3 

Summary  of  Mass-Sampling  Probe  Evaluation 
Tests  with' Refrigerant-11 


Test 

MacH 

No. 

Mixture  Total 
Temperature  (*R) 

Liquid  Input 
Mass  Percent 

Probe 

Location* 

Sampled 
Mass  Percen 

Error 

<!i) 

A 

<  1: 

750 

20.  6 

0.  33  R 

24.  8 

-J6 

13 

(pipe) 

29.  2 

0;  G7  R 

24.  2 

-17 

c 

62.  4 

0.  33  R 

57.  1 

-  9 

D 

G2.  1 

0.  67  R 

57.  4 

-  8 

E 

2.0 

COO 

30.8 

0 

42.5 

+  7 

P 

(nozzle) 

0.  58  R 

38- 6 

-  3 

G 

6C0 

•17.  2 

0 

4v.  8 

-  5 

II 

' 

0.  58  R 

45.2 

-  4 

I 

53.9 

0.  33  R 

50.  7 

-  6 

J 

2.0 

■100 

57.4 

0 

f 

75.0 

+31 

K 

(nc  heat) 

0.58  R 

65.6 

+  14 

I. 

67.  1 

0 

83.9 

+33 

M 

0.58  R 

78.5 

+  17 

Radius  fraction  from  centerline  in  pipe  or  nozzle  exit  plane. 
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When  the  mass -sampling  system  was  used  in  jet 
penetration  tests  in  which  unheated  a  ir-  was  used  (these 
results  are  described  later),  the  results  suggested  that 
the  sampled  mass  fraction  was  incorrectly  high.  To 
assess  these  findings,  evaluation  tests  J-M  in  Table  3 
were  conducted  with  unheated  air.  Indeed,  the  sampling 
errors  (+14  to  +33%)  were  greater  and  in  the  opposite 
direction  from  those  determined  with  heated  air.  In  addi¬ 
tion,  the  compositions  measured  at  0;  58  R  were  10  mass 
percent  lower  than  those  measured  at  the  centerline, 
which  suggested  that  the  test  stream  was  nonuniform,  due 
to  either  insufficient  mixing  or  condensation  on  the  pipe 
wall  surfaces.  If,  in  these  evaluation  tests  (and  the  jet 
^penetration  tests  described  later),  a  significant  portion  of 
the  injectant  remained  in  the  liquid  phase,  either  or  both 
of  two  effects^might.  cause  the -erroneously  high  results: 

(a)  liquid  drops  could  intermittently  plug  the  probe,  block¬ 
ing  the  air  flow  into  it;  and  (b)-if  the  shock  wave  on  the 
sampling  probe  were  detached,  the  centrifugal  forces 
created  by  the  shock  could  cause  the  lighter  air  component 
of  the  flow  to  be  deflected  by  a  larger  amount*  resulting  in 
disproportionately  large  amount  of  the  injectant  (both  liquid 
and  gaseous)  entering  the  probe. 

While  the  performance  of  this  probe  leaves  some¬ 
thing  to  be  desired,  the  qualitative  results  obtained  with 
it  are  of  considerable  interest,  as  we  shall  now  show. 

B.  DETERMINATION  OF  LIQUID  JET  FLOWFIELD 

Detailed  surveys  were  conducted  in  the  y-z  plane 
at  x  (s  x/d-)  =  50  for  refrigerant- 11  jets  that  were  injected 
normally  into  an  unheated  (T^  =  525°R),  Mach  2.  72  air- 

stream  with'  pro  =  1  atmosphere,  and  two  jet  conditions: 
d,  =  0.  027,  pf  ^  905;  and  d+  =  0.  052,  pf  =*•  302  (see  Table 

J  3  3 

2,  Tests  12  and  13).  Since  pressure  and  velocity 

gradients  exist  in  the  flow,  significant  differences  can 
exist  between  composition  distribution  and  the  distribu¬ 
tion  of  mass  flux,  which  is  a  property  of  main  interest. 
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Therefore,  pitot  and  cone- static  pressure  surveys  also 
were  made.  Pitot  and  cone-static  pressure  data  also 
were  measured  for  dj  =  0.  027,  pt.  =- .21 1  (Table  2,  Test  11), 
but  complete  mass-sampling  data ''were  not  obtained  for 
this  condition  because  it  was  established  that  a  sizable 
portion  of  the  injectant  occurred  below  y  =  4.  6  (y  =  0. 125 
inch),  which  was  the  minimum  value  at  which  data  could 
be  ^measured  because  of  the  size  of  the  probe.  Limited1 
pitot-pressure  data  were  measured  at  x  =  25  for  Test  13. 

Vertical  distributions  of  mass  composition  at  five 
lateral  distances  from  the  jet  centerline  are  shown  in 
Fig.  14.  Data  could  be  measured  only  down  to  y  =  4.  6 
and  2.  3  in  Tests  12  and  13,  respectively.  In  general,  the 
data  are  rather  well-behaved  in  regions  where  the  injec¬ 
tant  mass  content  is  less  than  ~  50%,  but  show  more  scat¬ 
ter  at  higher  compositions,  which  is  not  surprising  in 
view  of  the  probe  evaluation  test  results  (Table  3)  for 
comparable  conditions  of  stream  temperature  and  compo¬ 
sition  level.  Accordingly,  the  data  curves  drawn  for  the 
high  composition  regions  are  open  to  question.  This  fact 
will  be  examined  later  for  the  results  shown  in  Fig.  14a 
by  comparing  the  injectant  mass  flow  rates  computed  from 
both  the  solid  and  dashed  data  curves  drawn  for  z  =  0, 

2.  5,  and  5.  0. 

Some  typical  pitot  and  cone-static  pressure  data 
are  shown  in  Figs.  15  and  16  (the  complete  pressure  data 
for  Tests  11-13  are  given  in  the  Appendix).  These  data 
were  measured  with  the  probe  rakes  that  were  continuously 
traversed  in  the  lateral  (z)  direction  for  various  levels 
above  the  plate  surface  (y).  The  data  were  then  cross- 
plotted  to  obtain  pressure  distributions  for  various  lateral 
planes.  The  data  for  Test  13  were  taken  at  more  closely- 
spaced  vertical  positions  to  obtain  a  more  definitive  mea¬ 
sure  of  the  flow  structure  that  was  shown  to  be  heeded  on 
the  basis  of  results  obtained  :in  Tests  11  and  12.  Of 
particular  interest  are  the  pronounced  pressure  decreases 
that  occur  in  the  jet  region.  It  will  be  shown  that  this 
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TEST 


TEST  12 


feature  of  the  pitot-pressure  data,  in  particular,  is-useful 
for  determining  an  approximate  measure  of  the  jet  bound¬ 
ary  geometry.  The  pitot-pressure  maximum  that  occurs 
outside  the  jet  structure,  e.  g. ,  at  y  =  19.  2  for  x  =  25  in 
Test  13  (Fig.  17a),  is  caused  by  the  interaction  bow  shock, 
whose -position  corresponds  to  this  maximum  pressure. 

On  the  other  hand,  the  interaction  shock  had  already  inter¬ 
sected  the  air  jet  boundary  prior  to  x  =  50  in  Test  13,  and 
the  pitot-pressure  decrease  in  the  vicinity  of  z  =  25  (Fig. 
15a)  occurs  at  the  air  jet  boundary., 

A  difficult  problem  here  is  the  need  to  specify 
local  mixture  quality  and  the  local  liquid-phase  velocity, 
which  may  differ  significantly  from  that  of  the  gas  phase. 
There  also  is  a  question  of  the  extent  to  which' the  pres¬ 
sure  measurements  per  se  are  affected  by  the  liquid  phase. 
For  simplicity,,  the  present  data  are  reduced  on  the  as¬ 
sumption  that  the  gas  flow  properties  can  be  measured  in¬ 
dependently  of  the  presence  of  the  liquid  phase.  Then,  the 
gas  flow  properties  are  computed  for  either  of  two  assump¬ 
tions  on  mixture  quality:  the  injectant  is  all  gas  or  all 
liquid.  The  gas -phase  stagnation  temperature  is  taken  in 
both  cases  to  be  that  of  the  initial  air  plenum  temperature, 
which  was  only  slightly  lower  than  that  of  the  initial  injec¬ 
tant  plenum  temperature.  This  last  assumption  is  reason¬ 
able  even  for  the  case  where  the  injectant  is  all  vapo'r, 
because  itisexpected  that  the  energy  required  to  vaporize 
the  liquid  would  have  to  be  selr-supplied  and,  hence,  would 
decrease  the  overall  temperature  of  the  injectant. 

With  the  foregoing  assumptions,  the  gas  flow  prop¬ 
erties  are  computed  by  the  procedure  detailed  in  the 
Appendix.  The  gas-phase  static  pressure  and  Mach  num¬ 
ber  were  deduced  from  the  pitot-pressure  and  cone- 
static-pressure  data  by -combining,  the  Rayleigh  pitot 
formula  with  solutions  of  the  Taylor-Maccoll  differential 
equation  for  supersonic  flow  around  cones  using  the  appro¬ 
priate  value  of  the  specific  heat  ratio  y.  In  particular, 
y  =  1.  4  was  used  for  the  case  of  no  vaporization  (all  air), 
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and  ywas  based  on  the  local  mixture  composition  for  the 
case  of  complete  vaporization..  Similarly,  the  gas -phase 
molecular  weight  was  based  on  either  pure  air  or  the  local 
mixture  composition.  The  adiabatic  energy  equation,  the 
gas-phase  equation  of  state,  and  the  sound  speed  equation 
were  then  used  to  compute  the  local  gas-phase  mass  flux, 
from  which  the  local  injectant  mass  flux  was  deduced.  It 
is  important  to  note  that,  so  long  as  only  the  injectant 
mass  flux  is  of  interest,  there  is  no  need  to  specify  the 
liquid-phase  velocity  in  the  case  of  no  vaporization. 

The  pressure,  Mach-number,  and  mass-flux  dis¬ 
tributions  jomputed  ;for  the  case  of  complete  vaporization 
were  cross-plotted  to  deduce  the  contours  shown  in  Figs. 
17  and  18.  The  results  shown  for  the  regions  below  y  = 

4.  6  and  2.  3  in  Tests  12  and  13,  respectively,  are  esti¬ 
mates.  One  means  of  judging  the  validity  of  these  re¬ 
sults  is  to  compare  the  injectant  mass- flow  rates  (Wj's) 
determined  by  integrating  graphically  the  mass  flux  dis¬ 
tributions  with  the  metered  input  values  (Table  4).  When 
we  consider  the  dashed  curves  shown  for  Test  12  in  Fig. 
14a,  which  favor  the  lower  data  points,  the  error  in  the 
integrated  Wj  is  reduced  from  58  to  34%;  thus,  the  de¬ 
duced  flow  rate  is  very  sensitive  to  the  interpretation  of 
the  composition  data. 


Table  4 

Error  in  Integrated  Mass  Flux  Compared  to  Metered  w. 


Test 

12 

Test  13 

Measurement 

w.(lb/s) 

3 

Error 

\v.(lb/s) 

3 

Error 

Metered  Input 

0.  099 

-- 

0.  199 

— 

Sampling  (solid  curves,  Fig.  14) 

0.  156 

+58% 

0.255 

+28% 

Sampling  (dashed  curves,  Fig.  14a) 

0.  133 

+34% 

-- 

r  - 
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DIMENSIONLESS  HEIGHT  ABOVE  PLATE 


(a)  TEST  12 


DIMENSIONLESS  DISTANCE  FR 
PLATE  CENTERLINE,  z 


FIGURE  18 
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When,  for  Test  12,  the  injectant  is  assumed' to  be 
all  liquid,  the  deduced  flow  rate  is  about  twice  as  large 
(~0.  3  lb/ s).  Nevertheless,  some  portion  of  the  injectant 
probably  still  was  liquidat  x  =  50  in  these  tests,  because 
both  and  T^  were  below  refrigerant-llts  boiling  tem¬ 
perature  (5,35°R),  and  the  local  pressure  was  only  Slightly 
less  than  refrigerant- 1  l’s  vapor  pressure,  The  impor¬ 
tance  of  knowing  mixture  quality  is  evident 

Nevertheless,  we  have  assumed  that  the  mass  flux 
contours  based  on  the  all- vapor  assumption  (Pig.  18)  are 
qualitatively  useful  and  have  used  them  to  derive  spatial 
distributions  of  the  injectant.  Figure  19  shows  that  for 
Tests  12  and  13,  respectively,  66  and  70%  of  the  mass  is 
found  below  the  half-height  of  the  jet,  and  75  and  °5%  of  it 
lies  within  the  half-width-(i.  e.,.,  the  central  50%)  of  the 
jet,  These  results  are  very -significant  for  scramjet  en- 
gine  design.  In  addition,  for  Test  12,  the  mass  percent? 
age  curve  peaks  in  a  lateral  plane  «  3  jet  diameters  from 
the  centerline.  This  observation  is  consistent  with  the 
previously  described  photographic  data  for  jet  spreading 
that  showed  mass  apparently  being  stripped  from  the 
periphery  of  the  jet.  The  failure  to  find  an  off-axis  maxi¬ 
mum  for  the  dj  =  0..052  results  of  Test  13  can  be  ex¬ 
plained’ by  the  fact  that,  if  such  a  maximum  existed,  it 
could  have  been  removed  by  the  mass  diffusion  and  mixing 
processes,  which  would  have  been  operative _for  a  longer 
period  of  time  (by  about  a  factor  of  2,  since  x  was  twice 
as  great)  in  this  test. 

The  value  of  the  present  results  would  be  enhanced 
if  the  mass  distribution  curves  were  found  to  be  similar 
on  some  dimensionless  basis.  This  possibility  is  ex? 
amineci  in  Fig.  20,  where  the  local  jet  height  yg  and’ 
width  Zg  have  been  used  to  nondimensionalize  distance. 
Because  of  the  varying  degree  of  error  that  is  known  to 
exist  in  the  mass  flux  data,  it  is  not  expected  tuat  the 
distributions  from  the  two  tests  in  either  y  or  z  would 
reduce  to  a  single  curve.  Nevertheless,  they  are  fairly 
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DIMENSIONLESS  DISTANCE  FROM' PLATE  CENTERLINE,  z/zB 
OR  DIMENSIONLESS  HEIGHT  ABOVE  PLATE,  y/yB 


PERCENTAGE  OF  MASS  BELOW  y/yB, 

OR  PERCENTAGE  OF  MASS  WITHIN  ±  z/zB 

FIGURE  20  SIMILARITY  OF  DIMENSIONLESS  INJECTANT  MASS  DISTRIBUTIONS 
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similar,  even  though  the  x's  differed  by  a  factor  near  2. 
This  result  suggests  that,  in  the^absence  of  further  data, 
these  curves  can  be  used  to  estimate  mass  distributions 
for  other  conditions, _using  the  empirical  Eqs.  (1)  and  (10) 
to  determine  y  and  z  . 

In  this«  respect,  it  is  of  interest  to  examine  the  ex¬ 
tent  to  which  the  mass-sampling  results  of  Fig.  18  agree 
with  Eqs.  (1)  and  (10),  which  were  derived  from  tthe  photo¬ 
graphic  results  (Table  5).  It  is  recalled  that  Eq.  (10)  re¬ 
fers  to  the  maximum  width,  which  is  assumed  to  occur  by 
x,«  150,,  whereas  the  data  are  for  x  =  50.  However,  it  is 
estimated  from  the  photographic  data  (e.  g. ,  the  =  2.  72 
results  in  Fig.  10)  that  zB  increases  by  less  than  10%  in 
goi?g  from  x  =  50  to  x  =  150.  Referring  to  Table  5,  we 
see  that  the  empirical  equations  predict  yjg's  and  zg's 
that  are,  respectively,  16-19%  and  19-21%  lower  than 
those  (at  the  0%  level)  from  the  sampling  data.  It  must 
be  considered,  however,  that  injectarit  fluid  near  the  jet 
boundary  occurs  in  the  form  of  a  low- concentration  vapor 
and/or  extremely  small  droplets  and,  hence,  would  not 
be  detected  in  a  photograph.  Therefore,  it  is  reasonable 
to  use  a  finite  composition  level  in  making  our  comparison. 

Table  5 

Comparison  of  Jet  Boundary  Locations 


Test  12 

Test  13 

Basis 

,yB 

Z:B 

yB 

ZB 

Empirical  Formulas 

21.  f, 

17.  6+ 

11.8 

13.  0+ 

Mass  Fraction  Data  (0%  level) 

26.5  ' 

21.  6 

14.  1 

16.5 

Mass  Fraction  Data  (10%  level) 

21.  7 

17.  6 

12.5 

13.  1 

Mass  Distribution  (95%  of  mass) 

20.6  ; 

16.4 

10.9 

10.  8 

Pitot-Pressure  Data 

22.  8  , 

17.  6 

13.  1 

14.  9 

"I* 

Maximum  jet  width  —  assumed  to  occur  by  x  « 150 
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If  the  10-mass-percent  level  is  arbitrarily  chosen,  ;the 
yg's  and  Zg's  agree  well.  Another  means  of  evaluating 
the  usefulness  of  the  empirical  predictions  is  to  consider 
the  amount  of  mass  for  which  they  account.  If  we  deter¬ 
mine  from  the  dash-dot  curves  inFig.  19  the  y  and  z 
values  that  account  for  95%  of  the  injectant  mass,  we  find 
that  they  are  lower  than  the  empirically  predicted  yg  and 
Zg.  Hence,  we  conclude  that  the  empirical  predictions 
are  reasonably  accurate,  accounting  for~  96%  of  the  in¬ 
jectant  mass. 

Pitot-pressure  data  also  are  useful  for  a  rough 
determination  of  the  jet  boundary  geometry,  as  shown  in 
Table  "  which  includes  yg  and  zg  values  at  x  =  50  esti¬ 
mated  i  jet  boundary  contours  (solid  curves  in  Fig. 

21)  define  ^oci  of  points  at  which  the  knee  of  the  pitot- 
pressure  distribution  (in  the  y-  or  z-direction)  occurs; 
i.  e. ,  points  at  which  the  pitot-pressure,  as  it  increases 
above  its  minimum  value  (which  occurs^in  the  region  of 
maximum  injectant  composition),  reaches  the  nearly  con¬ 
stant  level  that  exists  between  the  jet  boundary  and  the  in¬ 
teraction  shock.  Referring  to  Fig.  15,  one  sees  that 
there  is  some  arbitrariness  in  defining  the  knee  of  a  pitot 
pressure  curve.  Nevertheless,  the  jet  boundaries  so  de¬ 
termined  are  in  better  agreement  with  the  mass-sampling 
data  (0%  level)  than  are  the  results  from  the  empirical 
equations,  as  shown  in  Table  5.  Figure  21  also  includes 
a  jet  boundary  contour  based  on  pitot-pressure  data  at 
x  =  25  for  Test  13  (dash-dot  curve),  and  it  is  noted  that 
the  pressure  data  show  a  16%  increase  in  yg  in  going 
from  x  =  25  to  x  =  50,  whereas  Eq.  (1)  predicts  a  12% 
increase. 
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DIMENSIONLESS  HEIGHT  ABOVE  PLATE,, 


DIMENSIONLESS  DISTANCE 
FROM  JET  CENTERLINE,  z 


FIGURE  21  APPROXIMATE  JET  BOUNDARIES  DEDUCED 
FROM  PITOT-PRESSURE  DATA 
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SECTION  V 


■  THEORETICAL  STUDIES  OF  JET  PENETRATION 


A.  GENERAL  FEATURES  OF  FLOW  MODEL 

The  usual  procedure  for  predicting  thejtrajectdry 
of  a  liquid  jet  element  of  infinitesimal  length  ds  and.  mass 
m  (Fig.  22)  is  to  apply  the  normal  component  of  the  mo¬ 
mentum  equation  to  obtain: 


F  =  -  mV.2/R  =  -  p.Ads  V.2  da/ds  ,  (11) 

n  3  3  3 

where  R  =  ds/dtt  is  the  local  radius  of  curvature  of  the 
trajectory,  and  pj  and  the  magnitude  of  Vj  are  assumed 
constant  at  their  initial  values.  The  advantage  of  using 
local  coordinates  n,  s  is  that  Eq.  (11)  holds  for  either 
constant  or  varying  m,;  when  the  velocity  of  the  element 
is  always  along  s.  Neglecting  the  relatively  small  gravity 
force,  the  force  Fn  can  be  represented  by: 


sin  gl  wds  , 


(12) 


where  CD  is  the  local  drag  coefficient  of  an  element  in  a 
crossflow  (a  =  90°).  Equations  (11)  and  (12)  are  combined 
to  give,  after_introducing  dimensionless  quantities  (A  s 
A/Ajj  q  2  q j / q^ ;  s  s,s/d w  =  w/d.), 

(w/A)  ds  =  -  (trq/2CD)  da/sin^a  .  (13) 


To  proceed,  it  is  necessary  to  specify  local  values  (as 
functions  of  s  or  a)  of  w,  A,  and  Cp,  which  are  all  func¬ 
tions  of  jet  distortion  and  disintegration  and  are  strongly 
interrelated.  When  this  is  done,  Eq.  (13)  and  the  rela¬ 
tions  dx  =  ds  cos  a  and  dy  =  ds  sin  a  are  used  to  determine 
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BOW  SHOCK 
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»URE  22  SCHEMATIC  REPRESENTATION  OF  LIQUID  INJECTION 
INTO  SUPERSONIC  STREAM 
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the  centerline  trajectory  in' the  x,  y  plane.  (In  general, 
it  is  necessary  first  to  determine  s  =  s(a)  before  these 
relations  can  be  applied.  > 


Although  the  cross-sectional  shape  of  the  liquid 
element  probably  is  noncircular,  Cq  is  often  taken  to 
correspond  to  the  drag  on  a  right  circular  cylinder  and 
is  then  based  oh  either  theory  or  experiment.  An  example 
of  the  former  approach  is  the  Newtonian  drag  theory,,  in 
which  local  pressure  on  the  forward  portion  of  the  cyl¬ 
inder  is  computed  by  assuming  that  particles  striking  the 
surface  lose  their  surface-normal  component  of  momen¬ 
tum.  The  pressure  distribution  around  the  body  is  then 
found  to  be: 


For 


-  ff/2  ^(ps  ff/2. 


2 

C  =  C  cos  <p  , 

P  % 


For  ff/2  <  <p<  -  ff/2. 


C 

P 


=  0 


(14) 


For  an  inclined  cylinder  (a  s*  90°),  only  the  component  of 
the  freestream  flow  that  is  normal  to  the  cylinder  axis 
contributes  to  the  drag,  in  which  case  Cp  =  (p  -  pe)/ 
qe  sin^».  A  second  expression  for  the  force  Fn  is  now 
determined  to  be: 


F  =  2  f  p  cos  (0  (w/ 2)  ckp  ds  , 
n  0 

(15) 

=  (2/3)  wds  C  q^  sin2a  . 

P0 

Combining  Eqs.  (12)  and  (15),  we  find: 

Cn  =  (2/3)  C  .  (16) 

D  po 


Here  Cp^  is  the  stagnation-point  (<p=  0)  pressure  coeffi¬ 
cient.  In  "regular"  Newtonian  theory,  Cp^  =  2,  but  a 
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constant  value  for  CD  is  incorrect  for  compressible  flows. 

*0 

Thus,  in  the  "modified-  Newtonian  theory,  Cn  is  based  on 

p0 

the  actual  stagnation  pressure  that  can  generally  be  cal¬ 
culated  as  a  function  of  M^..  For  subsonic -flow,  the  isen- 
tropic  flow  relationships  give: 

Cp^  =  12/yuJ)  J  [1  +  (y  -  1)  Me#2/2]y/(y_1)  -  1  J  .  (17) 

For  supersonic  flow,  the  normal- shock  relationships  give: 
cn  =  (2/yMx2)  i{(y+l)/(2yMo>2-y  +  l»1/(y-i)[(y  +  i)M  2 /  2]* 1  ^  .  A  , 


When  Eqs.  (;1 7)  and  ( 1 8)  are  applied-  to  an  inclined  cyl¬ 
inder,  Mco  should  be  replaced  by  Moo  sin  a,  so  that,  rigor¬ 
ously,  Cpn  =  Cpn  (Mooj  a)  should  be  used  in  integrating 
Eq.  (13). 

Variation  of  the  Newtonian  with' M*  is  compared 
with  experimental  data  for  a  solid,  normal  (a  =  90°),  con¬ 
stant-area  cylinder  (Ref.  15)  in  Fig.  23  (the  approxima¬ 
tion  shown  there  is  discussed  later).  The  difference  be¬ 
tween  data  and  theory  is  generally  sizable,  except  for  the 
hypersonic  regime  (M,*,^  5),  where  the  agreement  between 
data  and  the  modified  theory  is'  very  good.  The  deviation 
at  lower  Mm's  reflects  viscous  and  base-drag  effects  that 
occur  in  the  real  situation.  However,  none  of  the  curves 
necessarily  applies  to  a  liquid  jet,  since  jet  disintegra¬ 
tion  and  distortion  could  lead  to  viscous  and  base-drag 
effects  that  differ  significantly  from  those  occurring  on  a 
solid  body.  Moreover,  the  application  of  Cp  for  a  right 
cylinder  to  a  curvilinear  deflected  jet,  incorrectly  implies 
an  equality  between  bow  shock  and  jei  centerline  shapes. 
Hence,  it  is  concluded  that  no  really  satisfactory  basis 
exists  for  predicting  the  drag  on  a  liquid  jet. 
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The  situation -with  regard  to  predicting  jet  disinte¬ 
gration  and  distortion  is  . even  less  satisfactory.  One  of 
the  more  comprehensive  studies  in  this  area  was  made 
by  Clark  (Ref.  16),  who  determined_the  change  in  normal¬ 
ized  jet  cross-sectional  area,.  A.  s  A/Aj,  by  measuring 
the  electrical  resistance  variation  along  the  length  of  a 
water  jet  (V^  =  65  to  430  ft/ s)  in  a  nitrogen  crossflow 
(Voo  =  15  to  155  ft/s),  from  which  he  derived: 


For  c  <  1,  A  =  1  ,- 

o 

For  (21,  A  =  [0.  00392  (7.  68  c  -  1)  (c  -1)  +  1] 


(19) 


2  2 

where  c  =  [1  +  (Vj/V^,)  ]  s  /q  is  called  the  breakup  parame¬ 
ter.  The  occurrence  of  c  as  an-  independent  variable  re¬ 
sults  from  a  simple  theoretical  model  of  jet  breakup  at 
high  Weber  numbers  in  which  it  Is  assumed  that  jet  break¬ 
up  is  controlled  by  the  rate  at  which  the  jet  spreads  trans¬ 
versely  to  the  oncoming  stream.  The  spreading  is 
assumed  to  be  effected  only  by  the  dynamic  forces  imposed 
by  the  gas  stream  and  the  inertial  forces  within  the  jet. 

In  fact,  Clark  reasoned  that  jet  width  should  vary  as  fol¬ 
lows: 


w  =••  (1  +  2c)  .  (20) 

Clark's  data  arid  the  correlation  given  in  Eq.  (19)  are 
shown  in  Fig.  24.  While  considerable  scatter  (up  to  ±25%) 
exists  in  data,  their  uniqueness  makes  them  important 
and  useful. 


Before  proceeding,  it  is  important  to  note  that,  in 
formulating  the  flow  model  as  presented  here,  the  separa¬ 
tion  region  that  is  experimentally  observed  to  exist  in  the 
vicinity  of  the  orifice  exit  (see  Fig.  22)  is  neglected. 

Hence,  it  is  to  be  expected  that  analyses  based  on  the 
above  formulation  will  predict  an  incorrectly  low  penetra¬ 
tion  for  the  initial  portion  of  the  jet  trajectory;  independently 
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U<UP  AS  A  FUNCTION  OF  BREAKUP  PARAMETER 
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of  other  assumptions  made  in  a  particular  analysis.  In 
fact,  this  will’be  found  to  occur  for  those  analyses  dis¬ 
cussed  in  this  report. 

B.  REVIEW  OF  EXISTING  ANALYSES 


The  studies  of  Forde,  et'-ai.  (Ref.  6)  and  Catton, 
et  al.  (Ref.  8)  present  two  approaches  to  the  theoretical 
prediction  of  jet  penetration  that  are  fairly  representative 
of  the  scope  of  existing  analyses.  In  both  analyses,  cen¬ 
terline  trajectory  equations  are  derived  for  arbitrary  in¬ 
jection  angles  which,  when  specialized^  to  normal  injection, 
can  be  expressed  as  follows  (de  =  d..  is  assumed): 


Wo+kj 


In  K„  +  (x 


x0)/K2  + 


3 

|tK3 


+  (x  - 


x0)/K2 


where  Kj  2(1X4  -  K2)/K4,  K32  (xq  +  K4)/K4,  and  where 
K2,  K4,  x0,  and  y0 ,are  constants  that  depend  upon  the 
assumed  models  for  jet  cross-sectional  distortion,  disin¬ 
tegration,  and  drag  (xq  and  yg  actually  define  the  initial 
centerline  location;  however,  they  may  be  multivalued  for 
a  given  trajectory  if  the  integration  of  Eq.  (13)  is  done 
piecewise  because  of  changes  in  the  functional  form  of 
w,  A,  or  Cq).  To  assess  the  analyses  of  Refs.  6  and  8, 
we  shall  compare  their  predictions  of  jet  centerline  pene¬ 
tration  v/ith  the  empirical  correlation  of  jet  boundary  pene¬ 
tration  given  in  Eq.  (i).  (The  comparison  is  reasonable, 
since  in  the  theoretical  models,  the  difference  between 
centerline  and  boundary  locations  is  small  —  half  the  local 
-jet  thickness  in. the  n  -  direction. )  Observing  that  Eq.  (21) 
generally  does  not  predict  yc  «  qO*  5,  we  shall  compare 
the  equation  with  Eq.  (1)  for  a  range  of  q. 


In  Ref.  6,  jet  disintegration  and  cross-sectional 
deformation  are  neglected,  and  the  jet  is  assumed  to  be  a 
constant-diameter  cylindrical  body;  hence,  w/A  =  1,  Xq  = 
yQ  =  0,  and  K4  =  K2  =  Trq/2CD.  In  addition,  it  is  assumed 
that  CD  =  CD  (y ,  Mw)  3^  Cj^  (e» )#  where  ts  based  on  the 


(21) 
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Newtonian  theory.  Hence,  is  either  1-.  33,  for  regular 
Newtonian  theory  or  £  1.  23  -»  »)  for  . 'the  modified  New¬ 

tonian  theory  with  y  =  1.  4.  The  .prediction  for  CD  =  1.23 
is  compared  with  Eq.  (T)  for  q  =  3,  30,  and  100  in  Fig.  25, 
where  it  maybe  seen  that  for  q  a  3,  the  predicted  pene¬ 
tration  exceeds  the-data  for  x  >  35.  Increasing  Cq  to  1.  33 
does  not  noticeably  improve  the  agreement.  In -fact,  a 
Cd  as  high  as  2.13  (the  maximum  experimental  value 
shown  in  Fig.  23)  is  required/o  obtain  only  fair  agree¬ 
ment  over  the  range  of  q  tested. 


The  analysis  of  Ref.  8  incorporates  experimentally 
measured  effects  of  jet  distortion  and  disintegration.  They 
use  the  data  of  Clark  (Ref.,  16)  as  a  guide  in  deriving  an 
expression  for  local  jet' area  A;  in  particular,  they  assume 
that  A  =  1  up  to  a  point  Sq,  beyond  which  A  =  2.  25  q2s“^; 
hence,  sQ  =  (1.  5q)0*  5.  The  foregoing  expressions  are 
compared  to  Clark's  data,  n.  Fig.  24.  To  account  for  jet 
deformation,  w  is  assumed  to  iccch  a  maximum  value  of 
1.  6  at  Sq,  beyond  which  w  =  e  q^s“fe>,  where  e,  f,  and  g 
are  constants.  They  arbitrarily  assume  g  to  be  4,  which 
makes  w/A  independent  of  s  and  'simplifies  the  integra¬ 
tion  of  Eq.  (T3).  They  use  e  =  1,  and  f  =  1/25,.  to  give 
the  best  fit  between  their  predicted.and  measured  trajec¬ 
tories.  (Actually,  their  choices  for  various  constants 
appearing  in  the  analysis  result  in  some  inconsistencies; 
for  example,  the  values  chosen  for  e  and  f  do  not  give 
w  =  1.  6  at  s  =  Sq.  )  The  net  result  of  these  various  as¬ 
sumptions  is  to  divide  the  trajectory  into  two  regions,  in 
both  of  which  the  jet  is,  effectively,  a  constant-diameter 
cylinder.  In  the  first  region,  Xq  =  yg  =  0,  and  K4  =  K2  = 
irq/2C£),  as  in  Ref.  6,.  while  in  the  second"  region  (which 
starts’ at  the  ^following  values  of  Xq,  y^): 


IV  2,1  on  1  •  24.0,5  J. 

Xq  =  <K4  +  1.  63  q  )  -  iK4  , 


yQ  =  K4  in  I  1  +  x0/K4  +  [<Xq/K4)2  +  2xq/K4]0,  5  1 
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REDUCED  JET  PENETRATION.  (y/de)  (qj/qoo) 


FIGURE  25  COMPAP  SONS  OF  EXISTING  THEORETICAL  PREDICTIONS 
TO  PRESENT  DATA 
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and  K2  =  (1.  32ff/CD)  q  .  In  both  regions,  CD  is 
based  on  the  modified  Newtonian  theory  corresponding  to 
hypersonic  (M#  ■-*»)  flow.  In  particular  *  they  use  Cq  = 

(2/ 3)  (y  +  3 )/(y  +  1),  or  CD  =  1.  22  for  y  =  1.  4.  In  effect, 
this,  model  increases  (when  q  a  4.  9)  the  drag  in- the  second 
region  by  the  factor  0.  38  qO-62.  Figure  25  shows  that 
the  model  agrees  fairly  well  with  experiment  for  q  £  30, 
but  agreement  is  poor  for  q=“  3.  Observe  that,  this  model 
gives  a  smaller  dependency  of  the  reduced  penetration  on 
q  and  that  the  dependency  is  reversed  from  that  predicted 
by  the  analysis  of  Ref.  6.  For  the  same  level  of  CD 
(CD«  1.  2),  the  penetration  predicted  by  this  model  is 
considerably  improved  over  that  predicted  in  Ref.  6,  a 
fact  that  is  presumably  due  to  the  inclusion  of  jet  disinte¬ 
gration  effects. 

G.  PRESENT  ANALYSIS 

The  theoretical  studies  considered  in  the  previous 
section  have  shown  the  importance  of  including  jet  disin¬ 
tegration  effects  in  the  flow  modeling.  While  this,  was 
done  in  Ref.  8,  the  amount  of  empiricism  introduced  into 
that  analysis  limits  the  versatility  of  the  results  to  an  ex¬ 
tent  that  makes  it  difficult  to  classify  and  assess  the  ef¬ 
fects  of  yarious  assumptions.  Therefore,  the  analysis 
does  not  permit  optimum  use  of  experimental  data  to  sug¬ 
gest  further  improvements  in  the  theoretical  model.  More¬ 
over*  it  does  not  permit  computation  of  the  injectant  mass 
distribution.  Accordingly,  a  new  theoretical  analysis 
was  developed  for  predicting  the  centerline  trajectory  of 
a  normally  injected,  disintegrating  liquid  jet,  such  that 
the  effects  of  varying  jet  breakup  rate,  Mw,  and  CD  can 
be  assessed  by  comparing  the  predicted  injectant  mass 
distribution  with  experimental  data. 

In  the  present  analysis,  Clark's  (Ref.  16)  experi¬ 
mental  data  (Fig.  24)  were  fit  with  the  following  expres¬ 
sion  for  local  jet  area: 

For  Osc  sc  ,  A  =  1. 

0505  (22) 

For  €>€q»  A  =  exp  [-  0(e  ’  -cQ  *  )]  . 
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For  liquid  injection  into  supersonic' flow,  Vj/Va<z  1,  so 
that  [1  +(Vj/VC0)2]  «  1,  €  s^/q,  which  is  assumed  here¬ 

after.  Figure  24  shows  that  with  eg  =  2  and  0  =  2. 1,  ,Eq. 
(22)falls  within  the  region  of  certainty  of  the  data.  An 
advantage  of  the  form  of  Eq.  (22)  is  that,  for  conditions 
assumed  in ’this  analysis,  the  momentum  equation>(Eq. 
(13))  can  be  integrated  in:  closed  form.  Moreover,  Eq. 
(22)  does  imply  a  physically  reasonable  description  of  jet 
breakup  that  is  shown  by  its  differentiation: 


dA 

ds 


q  °* 


(23) 


i.  e. ,  following,  the  jet  path,  the  local' disintegration  rate 
is  proportional  to  local  jet  cross-sectional  area,  in¬ 
creases  with  the  square  root  of  the  freestream  dynamic 
pressure  and  decreases  with  the  square  root  of  the  jet 
dynamic  pressure.  However,  it  will  be  shown  later  that 


€  0  and  /3  may  vary  with  q,  so  that  the  exact  dependency  of 
the  breakup"  rate  on  q  may  be  different  from  q"0, 


Since  no  accurate  measure  of  the  jet  shape  is 
known,  we  make  the  simplifying  assumption,  as  others 
have  done,  that  the  jet  remains  circular  in  cross  section, 
so  that  w/A  =  A"0*  Finally,  we  assume  that  Cp  = 

CD  (y,  Ma)  t  C  £)(a),  .and  Cjj  is  given  by  the  modified  New¬ 
tonian  theory  with  y  =  1.4.  With  these  assumptions,  Eq. 
(13)  can  be  integrated  to  obtain  the  following  expressions 
for  jet  centerline  trajectory  for  normal  injection: 


For  90°  2  #  a  ocq,  x  =  Kg  (cosec  a  -  1) 

yc  =  Kg  cot(a/ 2)  .  v 

For  0  <  V  X  =  Kg(cosec  ap  -  1)  +  -  K7E(tt,a()) 


y  =  K  in  cot(an/2)  +  E(a,e:  )  , 

CO  u  u 
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For 


where 

30  -('o'O0'  5;  *0  =  COt’1<S0/K5) 

/  \ 

E(«,a0)  =  K6/K8  in  [{K?  +  Kg  -  tan  |)  (K?  -  Kg  -  tan  —)/ 

a. 

(K?  “  Kg  -tan|)  (K?  +Kg  -  tan  -jp)]; 

Kg  =  7rq/2GD;  Kg  =  (2/0). q°*  5;; 

K ?  =  'Kg  -  s0)/.K5;  Kg  =  (K?2  +  1)°*  5  . 


(Note  that  Eqs.  (24)  may  be  combined  to  give  Eq.  (21), 
with  x0  =  yQ  =  0. ) 

A  major  purpose  of  this  analysis  is  to  provide  a 
means  for  determining  the  downstream  mass  distribution 
of  the  injectant.  If  mass  torn  from  the  jet  is  assumed  to 
have  negligible  y-component  of  momentum,  so  that  fluid 
particles  remain. at  a.  nearly  constant  y  for  subsequent 
values  of  x,  the  y  distribution  of  mass  can  be  computed 
from  the  theory.  The  fraction  of  mass  G(a)  loss  at  a  (i.  e. , 
between  a  and  a  +  da)  is  given  by: 


For  90°  a  a  a  G(a)  =  0  .. 
a  <  a0,  G(a)  =  (2/Kg)[(K5/Kg)  cot  a  -TsQ/Kg)  +  1] 


(26) 


In  addition,  the  cumulative  mass  distribution  Ga,  where 
Gft  is  the  total  -mass  loss  between  ^  =  90°  and  a,  is  given 
by: 


For  90°  £  a  £  G  =  0  , 
u  a 

For  a  <  a _,  G  =  1  -  (K „/ 2)  G(a)  . 
0  a  6 


(27) 
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Both  Eqs.  (26)  and  (27)  can  be  combined  with  Eq.  (25)  to 
obtain  the  corresponding  distributions  in  terms  of  y,  but 
this  is  difficult  to  do  explicitly.  For  completeness,  we 
also  include  the  expression  for  s  =  s(a): 

For  90°  a-  a  s  =  Kc  cot  a  ,, 

°  5  (28) 

For  et  <aQ,  s  =  sQ  +  Kg  in  [1  +  (Kg/Kg)  (cot  a  -  cot  aQ)]  . 

Now,  since  Clark's  data  apply  to  liquid  injection 
into  subsonic  flow,  it  might  be  reasoned  that  q  in  the  fore¬ 
going  equations  should  be  more_properly  interpreted  as 
being  the  local  value  of  q'  h  q./qg  =  q/qg,  where  qs  is  the 
gas  stream  dynamic  pressure  behind  the  jet  interaction 
shock  wave,  so  that  the  dynamic  pressure  loss, across  the 
shock  would  be  accounted  for.  Tins  approach  causes 
some -difficulties,  however,  since  qg  depends  upon  the 
local  wave  angle,  which  decreases  from~  90°  (normal 
shock)  in  the  vicinity  of  the  injector  exit  down  to  the  Mach 
angle  corresponding  to  M*.  Furthermore,  the  deflection 
of  the  freestream  flow  by  the  wave  to  an  extent  that  also 
varies  with  the  local  wave  angle  means  that  qs  would  have 
to  be  further  modified  if  only  the  component  of  qg  in  the 
original  freestream  direction  were  of  interest.  For  sim¬ 
plicity,  however,  it  might  be  assumed  that  qg  and  hence, 
qg,  correspond  to -the  normal-shock  condition  for  Mw. 

Then  (Ref.  17), 

qs  =  [<y  -  l)/(y  +  D]  +  2 l(y  +  DM"  ,  (29) 

and,  with  the  exception  of  the  expression  for  K5,  q  in  the 
foregoing  trajectory  equations  is  replaced  by'(q/q  )  (the 
q  that  occurs  in  the  expression  for  K5  relates  to  the  drag 
coefficient,  which  already  accounts  for  the  fact  that  the 
flow  is  supersonic).  The  effect  of  this  assumption  on  the 
theoretical  results  is  examined  later. 

Before  proceeding  to  test  our  theoretical  model, 
it  is  of  interest  to  examine  the  effect  of  using  an  ar 


-  60  - 


dependent  drag  prediction.  This  can  easily  be  done  for  the 
case  of  a  constant-area,  nondisintegrating  jet  using  the 
following  expressions  for  the  drag  on  a  right  circular  cyl¬ 
inder: 

For  0  £  M  si.O,  CL  =  0.  667  (1  +  0,  250  M  2)  , 

co  oo  ' 

For  1.  0  <  M  ^  1.  46,  C^  =  0.  506  +  0.  328  M  ,  (30) 

oo  Ij  CO 

For  M  >  1.  46,  CL  =  1.  226  -  0.  51'4/M  2  . 

OO  jj  00 

It  is  shown  in  Fig.  2'3  that  the  above  equation  closely 
agrees  with  the  exact,  modified  Newtonian  drag  function. 
For  an  inclined,  cylindrical  jet  element,  M^  in  Eq.  (30) 
should  be  replaced  by  Ma  sin  a,  provided  it  is  assumed 
that  the  interaction  shock  shape  is  identical  to  the  jet 
boundary  shape.  The  .,  when  Eq.  (30)  is  substituted  into 
Eq.  (13),  with  w/A  =  1,  expressions  for  the  jet  centerline 
trajectory  can  be  determined.  These  expressions  are 
not  given  here  because  they  are  lengthy  and  change  their 
form  depending  upon  Mffl,  and  they  can  easily  be  generated 
with  the  use  of  available  integral  tables.  However,  Fig. 

26  compares  some  trajectory  calculations  based  on  them 
with  those  based; on  a,  drag  model  in  which  CD  *  CD  <MJ 
(i.  e. ,  no  a  dependency)  is  assumed.  It  is  seen  that  the 
effect  of  the  o-dependent  model  is  to  increase  the  pre¬ 
dicted  jet  penetration  at  the  low  Mach  number  (Ma  =  1'.  6) 
at  low  q  (q  =  3).  For  q=  30,  the  effect  is  smaller,  and 
for  q  =  100,  it  is  negligible.  For  the  higher  Mach  num¬ 
ber  (Mro  =5),  the  effect  is  negligible  even  at  low  q. 

D;  COMPARISON  OF  PRESENT  ANALYSIS  TO  DATA 

In  Fig.  27a  the  present  analysis  is  -compared  to 
the  empirical  Eq.  (1)  for  jet  boundary  penetration  (dashed 
curve)  using  e  Q  =  2  and  £  =  2.  1  (these  values  were  pre¬ 
viously  shown  to  give  a  good  fit  to  Clark's  data).  In  gen¬ 
eral,  the  predicted  curves  fall  far  below  the  curve  given 
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FIGURE  26  EFFECT  OF  DRAG  MODEL  ON  PREDICTED  JET  PENETRATION 
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by  Eq.  (1),  suggesting  that  the  disintegration  rates  used  are 
excessively  high,  with  ai  dependency  on  q  that  was  not  seen 
from  the  photographic  data,  and  it  is  noted  that  this  q- 
dependency  is  the  reverse  from  that  predicted  in  Ref.  8. 
Finally,  Fig.  27a,shows  a  slight  reduction  in  penetration 
when  Mm  is  increased  from  1‘.  6  (solid  curves)  to  5  (dash  - 
dot  curves).  This  result  is  interesting,  since  the  photo¬ 
graphic  results  shown  in  the  last  part  of  Fig.  6  indicate 
a  sizable  effect  in  this  direction-Creduced  penetration)  for 
a  smaller  change  in  Mffl  (increase  from  1.  62  to  2.  72).  In 
view  of  the  limited  accuracy  of  the  photographic  data, 
however,  it  is  felt  that  more  extensive  tests  are  required 
to  establish  firmly  the  effect. 

If  eg  is  kept  at  2,  but  /3  is  decreased  to  0.  65  (the 
variation  of  A  with  s  predicted  with  this  value,  is  com¬ 
pared  to  Clark's  data  in  Fig.  24),  an  improved  overall 
fit  between  theory  and  data  is  obtained  as.  seen  i'n  Fig.  27b. 
Although  the  discrepancy  at  low  q  is  reduced,  a  larger 
effect  of  q  is  predicted,  as  well  as  a  slightly  increased 
dependency  on  M^.  Figure  27c  presents  the  theoretical 
predictions  for  the  case  where  the  disintegration  rate  is 
based  on  conditions  existing  behind  the  interaction  shock, 
again  using  ,  =  2  and  j3  =  2.  1.  The  overall  spread  of 
the  various  curves  about  the  curve  from  Eq.  (1)  is  no 
better  than  in  part  (b),  and  the  effect  of  at  a  given  q 
is  now  reversed  (opposite  to  the  trend  seen  in  Fig.  6), 

From  these  results  we  conclude' that  and  /3 
must  somehow-'depend  on  q.  At  this  point,  one  could  per¬ 
form  a  numerical  experiment  to  ascertain  these  dependences 
on-q.  However,  since  we  are  more  concerned  with  predict¬ 
ing  the  cumulative  injectant  mass  distribution  Gy,  we  chose 
to  make  use  of  th  ;  available  experimental  Gy's. 

For  this  purpose,  (it  is  expedient  to  set  e  q,  whose 
effect  on  the  predicted  penetration  is  secondary  to  that  of 
0,  at  a  fixed  value  for  given  conditions  of  injection,  and 
then  seek  the  value  of  $  that  gives  the  best  agreement 


between  theory  and  data.  This  value  of  Cq  was  chosen  to 
reproduce  the  point  on  the  y-axis  intercepted  by  a  line 
that  (a)  has  the  same  slope  as  that  at  the  inflection  point 
(i.  e. ,  d^Gy/dy2  =  0)  of  the  experimental  Gy  curve  and  (b) 
passes  through  the  inflection  point.  The  values  of  y  de¬ 
termined  visually  in  this  manner  are  1.  3  for  q  =  3. 1 
(Test  13)  and  4.  2  for  q  =  10.  3  (Test  12)  (see  Fig.  28). 

The  corresponding  values  of  C  q  are  0.  55  and  1.  75,  re¬ 
spectively,  when  the  disintegration  model  is  based  on  q, 
and  0.  15  and  0.  50  when  the  model  is  based  on  q'. 

Using  the  above  values  of  Cq,  various  values  of 
P  were  tried.  Figure  28  shows  the  theoretical  Gy-distri- 
butions  for  narrow  ranges  of  p‘  to  show  their  sensitivity 
to/).  With  an  appropriate  /J,  agreement  is  good  in  the 
region  between  20  and  80%.  For  the  q  =  10.  3  data  in  Fig. 
28,  Cq  =  1.  75  and  £  =»  0.  30  give  a  good  fit  when  the  break¬ 
up  model*  is  based  on  q.  The  values  of  Cq  and  0  based  on 
Clark's  data  (2  and  2.  1)  give  a  poor  fit,  which  is  only 
slightly  improved  when  0  is  decreased  to  0.  65.  For  the 
q  =  3.  1  data,.  Cq  =  0.  55  and  0.  20  give  the  best  fit  for 
the  q-based  breakup  model.  The  vaiues<of  Cq  and  /)  de¬ 
termined  for  the  q  =  10.  3  data  do  not  give  as  good  a  fit 
here,  which  further  suggests  that  e0  and./)  are  functions 
of  q.  As  shown  in  parts  (c)  and  (d)  of  Fig.  .28,  the  break¬ 
up  model  based  on  q'  gives  about  equally  good-agreement 
when  the  previous  values  of  Cq  are  reduced  by  the  ratio, 
q/q',  and  different  values  of/J  are  used;  Again  the  agree¬ 
ment  is  poor  for  Cq  =  2  and  £  =  2.  1,  and  it  is  degraded 
when  an  attempt  is  made  to  fit  both  the  q  =  3.  1  and  10.  3 
data  with  the  same  values  of  c  q  and  /3. 

The  sets  of  C  q,  /3  that  hc.vc  been  shown  to  give 
good  data  fits  (1.  75,  0.  30  for  q  =  10.  3;  0:  55,  0.  20  for 
q  =  3, 1;  0.  49,  0.  55*  for  q'  =  36.  9;  and  0.  15,  0.  35  for 
q1  =  11.  1)  can  now  be  used  in  an  attempt  to  formulate  a 
dependency  of  c  q  and  0  on  q  or  q'.  Taken  at  face  value, 
the  above  results  suggest  that'CQ  «  q0.  96^  q0.  35^ 

c  q  «  (q')®’  an(j  p  «  (qi  )0.  45^  jiie  q  -  i o.  3  data  were 
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FIGURE  28  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL 
ACCUMULATIVE  DISTRIBUTIONS  FOR  INJECTANT  MASS 


used  to  determine  the  proportionality  constants  in  the  fore¬ 
going  expression^  When  this  was  done,  it  was  found,  by 
comparing'the  predicted  and  empirical  jet  trajectories, 
that  the  predicted,  values  for  Cq  were  excessively  high  for 
large  q  or  q'.  In  addition,  the  values  of  p  appeared  to  be 
too  low  for  large  q  or  q>,  particularly  for  the  q-based  0. 

In  view  of  these  results,  and  based  on  the  experience 
gained  in  judging  the  effects  of.  Cq  and  £  on  the  predicted 
‘trajectories,  the.  following  functional  fqrms  of  cQ  and  /3 
were  assumed:  *  /3  «  q°*  5  and  c  ^  pi  «  (q')°*  5.  To 

maintain  agreement  with  the  experimental  resultsTor  Gy, 
the  q  =  10.  3  data  were  used  to  compute  the  proportionality 
constants  which  resulted  in  the  following  equations: 

€Q  =  0.545  q°-5;i3  =  0.  094  q0,  5  , 
and  (31) 

€ 0  =  0.082  (q')0,5;  /3  =  0.091  (q»)°i5  . 

The  trajectories  predicted  with  the  use  of  Eqs. 

(31)  are  shown  in  Fig.  2^  for  q  values  between  3  and  100 
and  Ma,  =  1.6  and  5.  0.  The  overall  improvement  com¬ 
pared  to  Fig.  27  is  obvious  for  both  breakup  models.  The 
spread  of  the  theoretical  trajectories  about  the  empirical 
curve  (generally  high  for  q  s:  30,  low  for  q  =  3)  is  of  the 
same  order  of  magnitude  as  the  scatter  of  the  photo¬ 
graphic  data  on  which  the  empirical  equation  is  based 
(see  Fig.  7).  Hence,  there  is  no  firm  basis  for  much 
further  evaluation  of  the  theoretical  results  m  terms  of 
their'  predicted  dependency  on  q  and  Mw.  Moreover, 
there  is  little  to  choose  between  the  breakup  models. 

For  the  q-based  model,  the  predicted  penetration  is 
■nearly  independent  of  Mw  for  q  >  30,  whereas  for  the  q'- 
based  model,  it  increases  with  Mm  for  q  >  30.  We  do  not 
have  a  sufficient  basis  for  judging  whether  this  may  be  a 
•point  against  the  q1 -based  model. 

We  have  already  seen  in  Figs.  28a  and  c  the  ex¬ 
tent  to  which  the  values  of  Cq  and  0  given  by  Eqs.  (31) 
give  a  fit  to  the  q  =  10.  3  data.  The  extent  to  which  the 
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REDUCED  JET  PENETRATION,  yq 


DIMENSIONLESS  DISTANCE  FROM  INJECTOR,  x 

FIGURE  29  COMPARISON  OF  EMPIRICAL  AND  THEORETICAL  JET  TRAJECTORIES 
USING  IMPROVED  JET  BREAKUP  MODEL 
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q  =  3. 1  data  are  also  fitted  is  shown  in  Fig.  3Q.  The  pre¬ 
dicted  Gy  at  a  given  y  tends  tc^be  lower  then  the  experi¬ 
mental  value  except  near  the  jet  boundary.  However,  the 
fit  is  considered  reasonably  good,  especially  when  the 
experimental  difficulty  in  sampling  is  remem!  sred. 

Up  to  this  point,  we  have  considered  only  the  ex¬ 
perimental  distributions  for  G  ,  without  considering  the 
distributions  for  G(y),  the  percentage  of  mass  at  a  level 
y,  because  the  Gy  distributions  are  less  sensitive  to  the 
exact  details  of  tne  jet  breakup  model  and,  hence,  aie 
more  useful  for  comparing  with  experiment.  However, 
Fig.  3.1  compares  the  theoretical  distribution  for  G(y), 
computed  using  the  q-based  breakup  model  (the  result 
obtained  with  the  q* -based  model  is^  nearly  identical),,  with 
the  experimental  result  for  the  q  =  3. 1  case.  The  effect 
of  assuming  that  jet  bres;kup  is  delayed  for  a  finite  dis¬ 
tance  is  clearly  shown,  by  the  disagreement  at  small  y, 
where  the  greatest  difference  between  theory  and  experi¬ 
ment  exists. 


An  important  question  that  must  be  answered  now' 
is  wh  ,ther  the  jet  area  variations  predicted  by  .Eqs,  (31) 
...agree  with  Clark's  data.  Figure  32  shows  that  for  q1  z 
100  (or  q>  30),  the  q' -based  breakup  model  predicts 
values  for  A  that  are  consistent  with  the  data,  and  it  falls 
down  only  at  low  q,  whereas  the  q-based  model  gives 
poorer  agreement  at  all  q  levels.  Therefore,  the  q'- 
based  model,  which  approximately  accounts  for  the  dy¬ 
namic  pressure  losses  across  the  interaction  shock,  is 
judged  to  be  the  more  appropriate  one  *o  use,  especially 
if  subsonic  data  such  as  Clark's  are  to  be  ^elated  to 
supersonic  data.  The  ability  to  derive  a,  V  akup  model 
based  on  the  present  data  that  is  also  cons:  ,ent  with 
previous  data  i?s  particularly  satisfying,  since  this  gives 
con  idence  m,  at  least,  the  qualitative  value  of  the  pres¬ 
ent  data. 

We  have  now  established  that,  excepting  low 
values  of  q  (q^  3),  good  agreement  exists  between  theory 
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FIGURE  30  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL 
ACCUMULATIVE  DISTRIBUTION  FOR  INJECtlANT  MASS 
USING  IMPROVED  JET  BREAKUP  MODEL 
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DIMENSIONLESS  HEIGHT  ABOVE  PLATE, 


DIMENSIONLESS  JET  AREA, 


(and  data,  in  terms  of  either  predicted  jet  boundary  loca¬ 
tion  or  injectant  mass  y-distribution.  it  is  desirable 'to 
consider  possible  explanations  for  the  discrepancies  at 
low  q,  which,  primarily,  have  appeared<as  low  predicted 
boundary  locations.  One  feature  of  the  injection  process 
that  has  not  been  considered  heretofore  is  the  mass  diffu¬ 
sion  and  turbulent  mixing  effect  which  obviously  con¬ 
tributes  tJ  the  dispersion  of  the  injectant.  particularly  in 
the  jet  boundary  region.  Perhaps  this  effect  is  relatively 
more  important  for  a  low-q  jet,  which  has  lower  absolute 
values  of  penetration  than  a  higher-q<  jet. 

A  final  observation  concerning  the  jet  breakup 
model  deals  with  the  rate  expression  given  in  Eq.  (23): 

dA/ds  =  -/S  q  "*A  .  (23) 

With  /3  «  q0,  5,  as  assumed  herein,  the  above  expression 
reduces  to  dA/ds.  *  A,  i.  e. ,  following  the  path  of  a  jet 
element,  the  disintegration  rate  depends  only  on  the  local 
cross-sectional  area.  Intuitively,  it  seems  that  dynamic 
conditions  of  the  jet  and  the  airstream  should  affect  the 
disintegration  rate.  The  failure  of  the  finalized  form  of 
the  theory  to  predict  such  effects  is  taken  ns  another  indi¬ 
cation  that  the  q-dependency .assumed  -fpr  various  parame¬ 
ters  might  be  incorrect.  There  is  certainly  a  need  to  ob¬ 
tain  considerably  more  detailed  flowfield  data  for  a 
variety  of  conditions  to  establish  a  basis  for  improving 
the  theory. 
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SECTION  VI 

SUMMARY  AND  CONCLUSIONS 


The  dispersion  of  a  liquid  jet  ,that  was  injected 
transversely  from  a  flat- plate  surface  into  a  supersonic 
airstream  was  studied  by  photographic  and  in-stream- 
probe  measurements.  The  data  were  used  to  derive 
empirical  correlations  for  predicting  jet  boundary  pene¬ 
tration  and  spreading  and  as  the  basis  for  the  development 
of  a  theoretical  model  that  gives  a  satisfactory  predic¬ 
tion  of  the  transverse  distribution  of  jet  niass. 

The  photographic  measurements,  which  were 
based  on  back-scattered  light,  were  obtained  for  the  in¬ 
jection  of  water  and  isopropyl  alcohol  from  orifices  with 
diameters  of  0.  017,  0.  027,  and  0.  052-inch  into  air- 
streams  having  Mach  numbers  of  1.  62  and  2.  72.  Jet/ 
freestream  dynamic  pressure  ratio,  \q),  jet  Reynolds  and 
Weber  numbers,  and  freestream  Reynolds  number  were 
varied  by  factors  of  8,  ,10,  11,  and  3,-  respectively,  in 
the  jet  penetration  tests,  and  by  factors  of  23,  4,  9,  and 
5  in  the  jet  spreading  tests.  The  penetration  results  are 
correlated  to  within  ±  10%  by  an  empirical  equation  which 
predicts  that  dimensionless  penetration  (referenced  to  the 
effective  orifice' diameter)  at  a  dimensionless  downstream 
distance  varies  only  with  q-fJi  No  other  clearcut  depen¬ 
dency  of  penetration  on  either  the  Reynolds  or  Weber  num¬ 
bers  or  the  physical  properties  of  the  fluids  was  observed. 
The  dimensionless  maximum  jet  width  was  correlated  to 
within  ±  7.  5%  by  an  empirical  equation  which  predicts  a 
dependency  only  on,  q^*  25  (equally-good  agreement  was 
obtained  by  replacing  q  with  the  jet  plenum-to-free- 
stream  static  pressure  ratio).  The  empirical  correlation 
for  jet  penetration  developed  from  these  data  is  in  reason¬ 
ably  good  quantitative  agreement  only  with  existing  em¬ 
pirical  correlations  that  also  are  based  on  q;  the  agree¬ 
ment  is  within  the  ±  20%  variation  that  was  found  to  be 


possible  as  a  result  of  differences  in  photographic  tech¬ 
niques  that  could  occur  from  program  to  program.  Howr 
ever,  a  logarithmic  dependency  of  jet  height  on  down¬ 
stream  distance,  predicted  from  the  present  correlation, 
shows  improved  qualitative  agreement  with  the  data  as 
compared  to  the  previous  correlations,  which  generally 
predict  a  power-law  dependency.  The  jet  widths  deter¬ 
mined  from  the  present  correlation. are  only  |  to  |  as 
large  as  those  given  by  prior  correlations,  but  they  are 
substantiated  by  our  mass-sampling  data. 

For  the  purpose  of  sampling  a  heterogeneous 
stream,  a  resistance-heated  probe  and  a  heated  gas 
chromatograph  system  were  developed.  Although  evalua¬ 
tion  tests  of  the  probe  showed  that:  it  gave  an  incorrectly 
high  value  of  injectant  concentration  under  conditions  of 
highTiquid  loading  in  the  test  stream,  the  system  was 
deemed  to  yield  useful  qualitative  results.  Mass-sampling, 
pitot-pressure  and  cone-static-pi‘essure  surveys  were 
made  at  a  downstream  distance  of  50  jet  diameters  in  two 
tests  with  an  airstream  Mach  number  of  2.  72  and  with 
trichlorofluoromethane  as  the  injectant,  and  flowfield 
properties  were  deduced  on  the  assumption  that  the  injec¬ 
tant  occurred  as  all  vapor  at  the  sampling  location.  The 
total  injectant  flow  rates  deduced  from  the  data  exceed 
the  input  values  by  58%  in  one  test  and  28%  in  the  other; 
these  values  are  twice  as  high  when  the  injectant  is 
assumed  to  occur  in, the  liquid  phase.  Accepting  the 
qualitative  value  of  the  all-vapor  results,  they  show  that 
65  to  70%  of  the  mass  is  below  the  jet  half-height,  and 
75  to  85%  of  it  is  within  the  central  jet  half-width.  In 
addition,  the  cumulative  mass  distributions  for  either 
the  transverse  or  lateral  directions  from  the  two  tests 
are  nearly  similar  when  the  respective  distances  are 
referenced  to  either  jet  height  or  jet  width,  even  though 
the  absolute-downstream  distances  differed  by  a  factor 
of  2  (different  orifice  sizes  and  injection  pressures  were 
used  in  the  two  tests).  This  fact  suggests  that  the  pres¬ 
ent  data  are  useful  for  estimating  injectant  mass  distri¬ 
butions  for  conditions  quite  different  from  those  used  in 


the  tests.  The  sampling  data,  which  give  rather  accurate 
measures  of  jet  height  and  width,,  showed  that  the  p'.ioto- 
graphically-based  empirical  correlations,  summarized 
above  give  slightly  low  values,  but  that  they  account  for 
at  least  96%  of  the  .injectant  mass,  and' that  the  injectant 
that  occurs  outside  the  empirically  predicted  boundaries 
is  at  a  concentration  that  is  probably  too  low  to  be  de¬ 
tected  by  ordinary  photographic  -procedures. 

Existing  theoretical  treatments  that  include  the 
effect  of  jet  disintegration  give  reasonable  predictions 
for  jet  penetration,  but  they  neither  allow  amassessment 
of  the  effects  of  various  necessary  assumptions  nor  pro¬ 
vide ‘for  deducing  mass  distribution.  An  analysis  was 
developed  on  the  assumptions  that:  (a)  the  jet  is  a  circu¬ 
lar,  disintegrating  body,  (b)  its  drag  can  be  based  on  the 
modified  Newtonian  theory,  (c)  its  local  cross-sectional 
area  can  be  based  on  Clark's  (Ref.  16)  data  for  jet  disin¬ 
tegration  in  a  subsonic  stream,  and  (d)  particles  torn 
from  the  jet  remain  at  a  constant  transverse  position  for 
subsequent  downstream  locations.  The  predicted  trans¬ 
verse  distributions  of  mass  agree  reasonably  well  with 
the  data  when,  the  jet  disintegration  parameters  appearing 
in  the  analysis  \re  based  on  approximate  conditions  exist¬ 
ing  behind  the  j  l  interaction  shock.  The  predicted  jet 
boundary  trajectories  show  dependencies  on  q  and  free- 
stream.  Mach  number  that  were  not  established  from  the 
photographic  data,  but  these  effects  are  comparable  in 
magnitude  to  the  scatter  of  the  photographic  data. 

The  detailed  flowfield  data  have  proven  to  be  ex¬ 
tremely  valuable  for  describing  liquid  jet  disintegration. 
Additional  effort  is  required  to  improve  the  accuracy  of 
the  mass-sampling  probe,  and  additional  mass-sampling 
data  for  a  wide  variety  of  test  conditions„are  needed  for 
further  assessment  and  improvement  of  the  theory. 
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APPENDIX 


EXPERIMENTAL  PRESSURE  DATA  AND  FLOW FIELD 
DATA  REDUCTION  TECHNIQUE 


Figure  A-l  presents  the  pitot-pressure  data  for 
Tests  11-13  at  x/dj  =  50,  and  Fig.  A -2  presents  the  cor¬ 
responding  cone-static  pressure  data. 

FLOWFIELD  DATA  REDUCTION  TECHNIQUE 

The  method  used  to  deduce  flow-field  properties 
from  the  instream  probe  data  is  based  on  the  following 
assumptions: 

1.  The  experimental  probes  respond  only  to  the 
axial  motion  of  the  flow,  or,  equivalently,  the  motion  is 
one -dimensional, 

2.  The  gas  phase  flow  can  be  measured  and  com¬ 
puted  independent? v  of  the  presence  of  the  liquid  phase, 

3.  The  portion  of  the  injectant  that  vaporizes  is 
m  thermal  and  dynamic  equilibrium  with  the  air,  and 

*1.  The  mean  flow  is  steady. 

The  pressure  coefficient  of  the  cone-static-prcs- 
sure  probe  is  given  by: 

Cp  =  (pc  -  p)/(y PM2/2)  ,  (A-l) 

where  pc  is  the  average  of  the  four  cone  surface-pressure 
measurements.  Values  of  Cp,  which  are  determined  from 
solutions  of  the  Tavlor-Maccoll  differential  equation  for 
subsonic  flow  around  cones,  arc  given  in  Ref.  18  for  vari¬ 
ous  values  of  the  cone  half-angle  and  M,  but  only  for 
values  of  y  =  1.4  and  1.33.  Fortunately,  however,  Cp  is 
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DIMENSIONLESS  DISTANCE  FROM  JET  CENTERLINE 


FIGURE  A— 1  PITOT-PRESSURE  DATA 
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DIMENSIONLESS  DISTANCE  FROM  JET  CENTERLINE 


DIMENSIONLESS  DISTANCE  FROM  JET  CENTERLINE. 


(c)  TEST  13 
x  =  50 


0.5  0.3 


PITOT-PRESSURE/AIR  PLENUM  PRESSURE,  pt7p(00 


FIGURE  A-1  (CONCLUDED) 


quite' insensitive  to  y  (see  Ref.  19)and  the  values  applica¬ 
ble  toy  =  1.  4  are  used  herein.  The  measured  pitot-pres¬ 
sure  is  introduced  into  Eq.  (A-l),  and  the  equation  is  re¬ 
written  as: 

Pc/Pt*  =  <P./Pt* )[1  +  (VM2Cp/2)]  .  (A-2) 

A  relation  for;p/p^  is  given  by  the  Rayleigh-  pitot  formula 
(Ref.  17): 


2yM 2/(y  +  1)  -  <y-  l)/(y  +  l)\  1/(y_1) 


Uv  +  DM 


(A-3) 


With  Eqs.  (A-2)  and  (A-3)  and  the  values  of  C  as  a  func¬ 
tion  of  M,  a  plot  of  pc/pt'  as  a  function  of  y  and  M  can  be 
derived,  from  which  M  can  be  determined  once  y  is  spe¬ 
cified.  Then,  p  is  determined  from  a -plot  of  Eq.  (A -1 ), 
written  as  follows: 

p/po  =  [1 +(yM2C  /2)]”-  .  (A-4) 

c  p 


The  procedure  beyond' this  point  depends  upon  whether 
the  phase  of  the  injectant  is  assumed  to  be  all  liquid  or 
all  vapor.  (Actually,  one  could  consider  the more  gen¬ 
eral  case  where  only, a  portion  of  the  injectant  Was  vapor¬ 
ized;  then,  however,  it  is  necessary  to  specify  the  ve¬ 
locity. of  the  condensed  phase,  which  is  generally  unknown. ) 
We  shall  consider  these  two  cases  separately. 


ALL  VAPOR 


When  the  injectant  occurs  as  ail  vapor,  the  mass¬ 
sampling  measurement,  which  gives  the  local  mole  frac¬ 
tions  of  refrigerant-11,  X^,  and  air,  X^  -  1-X^,  can 
be  used,  to  determine  several  required  properties  of  the 
gaseous  mixture.  The  specific  heat  ratio  of  the  mixture 
is  given  by: 
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where  yA  =  1. 4  and  y^  =  1.136.  'The  molecular  weight  of 
the  mixmre  is: 

'■¥4^'  <A-6) 

wherem^  =  28.-965  andlTl^  =  137.38.  We  now  determine 
the  mixture  gas  constant,  R,  and  the  refrigerant- 11  mass 
fraction,  Y^,  to  be: 

R  =  1545/m  YR=XRrHRM.  (cA'r-7)- 

The  static  temperature,  T,  is  determined  from  the  energy 
equation  for  adiabatic,  flow  in  the  form: 

T/Tt  =  [1  +  (y  -  V)M2/2f 1  ,  (A -8) 


where  Tt,  the  mixture  stagnation  temperature,  is  assumed 
equal  to  the  initial  air  plenum  Temperature.  The  remain¬ 
ing  properties  of  the  mixture  are  calculated' from: 


Speed  of  sound  definition 

V  =  M(yRT)0,  5 

(A-9) 

Equation  of  state: 

P  =  p/HT  . 

(A -10) 

The  mass  fraction  of  the  injectant 

is  then  given  by: 

<pV)R  =  YR(pV>  . 

(A-ll) 

ALL  LIQUID 

The  injectant,  when  it  is  assumed  to  occur  as  all 
liquid,  is  further  assumed  to  have  a  negligible  influence 
on  the  dynamic  and  thermal  properties  of  the  gas  .phase 
(pure  air).  Then: 


y  =  =mA;  k  =  (A.I2) 

(  I 

,  _  / 

and  T-^  V,  and  pjare  calculated  as  before  to  determine  the 
air  mass  flux,  (f)V)^>  The  mass  fraction  of  refrigerant-1] 
Yr,  is  still  given  by  Eq.  (A-7)  w*  thirl  given  by  Eq.  (A -6)  " 

Thftn*  •  ' 


|h„  s 
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